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SUMMARY STATEMEN T OF CONTRACT RESEARCH. CONCLUSIONS. AND 
ACCOMPLI SHMENTS 

The contract research has been conducted in the following 
three major areas: 

1. Analysis of numerical simulations and parallel 
observations of atmospheric blocking. 

2. Diagnosis of the lower boundary heating and the 
response of the atmospheric circulation. 

3. Comprehensive assessment of long-range forecasting 
with numerical and regression methods. 

The essential scientific and developmental purpose of this 
contract research is to extend our capability of numerical 
weather forecasting by the comprehensive general circulation 
model. The systematic work as listed above is thus geared to 
developing a technological basis for future NASA long-range 
forecasting. 

Experiments have been conducted for a 46— day summer period 
from July 1 to August 15, 1979, using the high resolution 2 x 
2.5 degree version of the Goddard Laboratory for Atmospheres 
(GLA) general circulation model (GCM) , to simulate the summer 
blocking observed during the period in the Northern Hemisphere. 
During the simulation, sea surface temperature (SST) anomalies 
were updated daily with observations to provide a realistic 
surface heating field. A parallel control run used climatology 
SSTs, in lieu of observed SSTs. Both the SST updating 
experiment and the control run were able to simulate the 
persistent blocking patterns in Eurasia and the eastern Pacific 
during July. However, only the SST updating experiment was 
able to generate the realistic blocking episodes in the last 
10-day period during the 46-day simulation. Forecast skill was 
examined in terms of 500 mb anomaly correlation and error 
kinetic energy. When the blocking patterns develop or become 
reinforced, a high forecast skill is apparent in nearly half of 
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the Northern Hemisphere, including the area ot blocking and its 
downstream region. Unlike the major winter blockings, which 
may be treated as the manifestation of wave number one and two, 
the summer blockings are localized, and have a smaller scale 
than that of winter blockings. However, summer blockings 
treated in this study are persistent, and possess a 
considerable predictability. The effort in this simulation is 
detailed in Kung, Min, Susskind and Park (1992). 

in an attempt to examine the association of SST anomalies 
with major winter blocking, twelve major blocking episodes in 
the Northern Hemisphere winter are selected for examination 
from a 34-year data period: four episodes for each category of 

Pacific, Atlantic and double blocking. It is shown in Kung, 
Susskind and DaCamara (1993) that the major single blocking in 
the Pacific or Atlantic is formed through constructive 
interference of the traveling zonal wave of n=l and the 
stationary wave of n=2 . The n-1 is supported by the energy ^ 
input through the nonlinear wave-wave interaction, and the n 2 
by in situ warming over the Pacific and Atlantic. The 
concurrent double blocking in the Pacific and Atlantic is 
formed by the stationary n=2 when the traveling n=l is weak. 

The cooling of the Pacific and Atlantic is associated with the 
single blocking in the Pacific and Atlantic. On the contrary, 
the warming in both the Pacific and Atlantic is associated with 
concurrent Pacific and Atlantic blockings, indicating the 
baroclinic nature of the double blocking. The hemispherical 
and regional energetics of blocking areas are examined in terms 
of kinetic energy and nonlinear kinetic energy transfer in n 1, 

2 , and 3 . 

Large-scale modes of variations in SSTs and the 
tropospheric circulation are examined with major principal 
components utilizing monthly mean fields of the global SSTs and 
Northern Hemisphere geopotential height (Z) at 700, 500 and 300 

mb levels in Kung, Chern and Susskind (1994). The data period 
covered is from 1955 to 1992. It is found that the 
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heterogeneity of SST data due to availability of satellite 
observations and difference of analysis schemes may result in a 
large systematic bias in the dataset. However, the bias may be 
effectively corrected through elimination of an appropriate 
principal component, n=l in this case. The first three 
components of monthly SST and Z fields during the 38-year 
period are scrutinized. The El Nino-Southern Oscillation 
(ENSO) mode of variation is observed in both 1st and 2nd 
components of SSTs. The inter-annual variations of principal 
components of monthly SST and Z fields are utilized to probe 
the association of SST components and Z components. Further, 
cross-correlation patterns of principal components of SSTs in 
reference to the tropospheric circulation are studied with 500 
and 300 mb fields. It is shown that the tropospheric response 
to SST anomalies are consistent at 500 and 300 mb levels, and 
the seasonal— range predictability of the tropospheric 
circulation is recognized with SST anomaly fields. 

Long-range predictability of the Northern Hemisphere 
circulation with numerical-dynamical and empirical-regression 
®®thods has been critically evaluated in reference to prominent 
circulation phenomena, including blocking, monsoon, and others 
(Kung, Chern, Susskind and Lin 1992). Whereas the 
deterministic model forecasting can be accomplished to 1.5 
months with the current GCM, it appears necessary and desirable 
to formulate the multiple regression forecast scheme at the 
range of l month to 1 year. 

Accomplishments of this contract research are detailed in 
4. Technical Details of Accomplished Results. 

2. RECOMMENDATIONS 

It has long been stipulated that SSTs may affect the 
formation of blocking. It has been shown in this contract 
research that SSTs are an essential forcing function for 
blocking development. Blocking is the persistent feature in 
the global circulation which can be tagged as the basic 
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circulation pattern to be forecasted in the long-range. Other 
features of the general circulation may be considered in 
reference to the predicted blocking occurrence. 

This contract research specifically revealed that there is 
a significant predictability of blocking if the SST component 
of the global circulation is properly prescribed. In the model 
forecast, it should mean the development of the joint ocean- 
atmosphere model. We have indicated that at least a 1.5 month, 
and possibly longer, forecast is feasible with GCMs. For the 
seasonal-range forecast up to 1 year, which is not in the 
immediate reach of the GCM forecast, a development of the 
empirical regression scheme should be considered. The results 
obtained in this study offer useful information for further 
diagnosis of the SST and circulation field, which will be vital 
in the joint development of ocean-atmosphere model forecasting 
and empirical-regression forecasting. 
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4 . TECHNICAL DETAILS OF ACCOMPLISHED RESULTS 

Technical details of accomplished results are presented in 
the following form of publication reprints for 4.1, 4.2, and 
4.3, and in prepublication form for 4.4. 

4.1 An Analysis of Simulated Summer Blocking Episodes, page 6. 

4.2 Long-Range Forecasting of the Northern Hemisphere 
Circulation with Numerical-Dynamical and Empirical- 
Regression Models, page 19. 

4 . 3 Prominent Northern Hemisphere Winter Blocking Episodes and 
Associated Anomaly Fields of Sea Surface Temperatures, 
page 47. 

4.4 Large-Scale Modes of Variations in the Global Sea Surface 
Temperatures and Northern Hemisphere Tropospheric 
Circulation, page 67. 
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4,1 An analysis of simulated summer blocking episodes 

By ERNEST C. RUNG and WEI MIN 

Department of Atmospheric Science, University of Missouri-Columbia, Columbia, Missouri 65211 

and 

JOEL SUSSKIND and CHUNG-KYU PARK* 

Laboratory for Atmospheres, Goddard Space Flight Center/NASA, Greenbelt, Maryland 20771 

(Received 8 April 1991; revised 30 July 1991) 


Summary 

Using the high-resolution 2 x 2.5 degree version of the Goddard Laboratory for Atmospheres general 
circulation model, experiments were performed for a 46-day summer period from July to 15 August 979. 
During the simulation, sea-surface-temperature (s.s.t.) anomalies were updated daily with observed s.s.t.s to 
provide a realistic surface heating field. A parallel control run used climatology s.s.t.s. in lieu of observed 

S S ' Both the s.s.t. updating experiment and control run were able to simulate the persistent blocking patterns 
in Eurasia and the eastern Pacific during July. However, only the s.s.t. updating experiment was able to 
generate the realistic blocking episodes in the last 10-day period during the 46-day simulation. Forecast skill 
was examined in terms of 500 mb anomaly correlation and error kinetic energy When the blocking patterns 
develop or become reinforced, a high forecast skill is apparent in nearly half of the northern hemisphere 
including the area of blocking and its downstream region. Unlike the major winter blockings, which may be 
treated as the manifestation of wave number one and two, the summer blockings are localized, and have a 
smaller scale than that of winter blockings. However, summer blockings treated in this study are persistent, 
and possess a considerable predictability. 


1. Introduction 

In a preceding paper Kung et al. (1990) reported that simulation experiments, using 
the high-resolution Goddard Laboratory for Atmospheres (GLA) general circulation 
model (GCM) with a realistic sea-surface-temperature (s.s.t.) field, were able to generate 
two successive major blocking events over the northern hemisphere Pacific and Atlantic 
oceans during January 1979. It was suggested that the development of the second blocking 
episode over the Atlantic toward the end of the month was due to the realistic heating 
field that provided an adequate baroclinic energy source for wave-wave kinetic-energy 
transfer, whereas adequate baroclinic energy existed within the initial data for the 
development of the Pacific blocking immediately following the initialization of the GCM. 
The formation of realistic blocking during the month was associated with the improved 
forecast skill of the ultralong waves. 

In the simulation experiments the ocean surface heating field was provided with 
daily updating of observed s.s.t.s (the procedure was developed by the third and fourth 
authors of this paper). One view on this series of experiments is that what was being 
tested is the upper limit of what would be possible with a coupled ocean-atmosphere 
model, since updating the observed s.s.t.s during the atmospheric-model integration may 
imply a correct ocean-model integration. Despite the encouraging results, however, a 
number of important questions were left unanswered in the preceding study. Only one 
initial state was involved in the experiments, and as such, the results cannot be gener- 
alized. The effects of uncertainties in initial conditions and the dependence of simulations 
on the initial regime of flow must be considered (see Palmer 1988, Tracton ct al. 1989, 
Tribbia and Baumhefner 1988). Although the improved blocking simulation resulted 
from the use of a realistic s.s.t. field, it was not determined whether this was due to the 

* Affiliated with General Sciences Corporation, Laurel, Maryland. 
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improved initial s.s.t. field or the real.stic time variations of s.s.t.s during the simulation. 
I he skill of blocking simulations in the northern hemisphere winter is often a major 
concern in medium-range forecasting, and it is of interest to investigate the simulation 
ot summer blocking episodes in similar GCM integrations. 

In the subsequent study two simulations during the summer of 1979 were performed 
Using an identical h.gh-resolution version of the GLA GCM with the GLA initial data 
f ' i° r '»». ■ model , meg, alien was carried out for ,he period Tm” 

15 August 1979. The s.s.t. anomalies were updated daily with observations in the same 
manner as in the preceding winter simulations. A parallel control run was performed in 
the same way, except that the s.s.t. climatology for the simulation period was utilized in 
ieu of observed s.s.t.s. Both the simulation with updated s.s.t.s and the control run are 
compared with observations. Blocking episodes and their patterns in the northern 
hemisphere are identified. The forecast skill is examined in terms of anomaly correlation 
and error kinetic energy at the 500 mb level. Considerable forecast skill is achieved in 
the summer simulation in association with blocking development and maintenance 

Effects of s.s.t. updating on the simulation are studied with the time variation of s s t 
anomalies during the period. 


2- Simulation experiments 

rr 46 ; da y simulations were performed using the 2 x 2.5 degree version of the 
GLA GCM for the period 1 July-15 August 1979. The GCM was the same high-resolution 
version as employed in Kung « al. (1990). The GLA fourth-order global atmospheric 
model is described in Kalnay-Rivas et al. (1977) and Kalnay et al. (1983). The hieh- 
resolution version is identical to the original coarse-resolution version with respect to 
physical parametrizations. The model was initialized with the GLA gridded analyses of 

Research pr<,8ram) oi ” w 

A blended global monthly s.s.t. analysis at the GLA based on in situ (ship and 
buoy) and satellite-retrieved data by the third and fourth authors of this paper was 

Router ^ in SuSSkind et aL < 1984 ) and Susskind and 

Reuter (1985). For the GCM experiment, daily gridded s.s.t. fields were obtained by 

linear interpolation of the blended monthly fields, assuming that the observed monthly 
means are located in the middle of each month. The interpolation in this manner is valid 
because of the smallness and slowness of the time variation in s.s.t.s. The GCM simulation 

l gl °f b H al ?TM Phere ’ T ith ^ S S t S updated durin8 the model integration using 
h.s set of daily fie ds, is referred to as the ‘experiment’ in this paper to contrast with 

the parallel control run and observed atmosphere. Daily s.s.t. fields obtained by the 
interpolation of monthly climatological s.s.t.s (Reynolds 1982) were utilized in the control 
run of the simulation. The Geophysical Fluid Dynamics Laboratory analyses of the 
FGGE observations during the simulation period were used to describe the observed 
circulation. The simulation and observation fields analysed in this study include twice- 
daily values of the circulation fields at 0000 and 1200 GMT. 


3. 


Observed sea-surface-temperature anomalies and blocking episodes 


p- ° bserved anomaly patterns of blended s.s.t.s from 30°S to 60°N are illustrated in 
g. 1 with four selected 5-day period averages during the simulation. The National 
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obtained by linear interpolation of the monthly fields, an inherent smoothing of the fields 
may obscure time variations of s.s.t.s. However, it may be suggested that the anomaly 
pattern during July is stable, and its time change very gradual. The pattern is characterized 
in the northern hemisphere by a large area of negative anomaly in the mid Pacific with 
positive anomalies to the north and south. The negative anomaly area in the mid Atlantic 
is small, and positive areas to the north and south dominate the Atlantic. Despite the 
possible oversmoothing of the s.s.t. variations, there are noticeable changes from the 
July to August pattern. In August the large negative area in the mid Pacific becomes 
small and fragmental, and the positive areas are enlarged and strengthened. In the 
Atlantic the negative and positive areas are both weakened in August. 

Distinct blocking episodes during the simulation period are identified in Eurasia and 
the eastern Pacific with the observed northern hemisphere circulation. They are listed in 
Table 1 for the identified periods, the subperiods of reinforcement, and the primary 
longitudinal locations in the 46— 66°N belt. The observed blocking episodes are compared 


TABLE 1 . Observed blocking episodes during the simulation period 
1 July-15 August 1979 


Period 

Subperiod 

Primary location 
(46-66°N) 

Eurasia: 



1-6 July 
6-26 July 

6-10 July 
11-15 July 
16-20 July 
21-26 July 

(90-120°E) 

(lO-SHTE) 

(10-90°E) 

(30-100°E) 

(75-140°E) 

31 July-4 August 
8-13 August 
10-15 August 


(90-130°E) 

(0-40°E) 

(75-120°E) 

Eastern Pacific: 



3-6 July 
13-17 July 
24 July-15 August 

24-28 July 
29 July-1 August 
2-4 August 
5-9 August 
10-14 August 

(150-90^) 
(140-100°W) 
(150-1 10°W) 
(150-1 lO'W) 
(150-110°W) 
(150-120°W) 


with the blocking occurrences identified in the experiment and control run in the 
longitude-time diagrams of Fig. 2. Because of the localized pattern of summer blocking, 
identification of blocking is done manually through observation of daily 500 mb geo- 
potential fields (see Kung and Baker 1986) rather than using an objective scheme such 
as performed by Kung etal. (1990). Figure 3 further illustrates the observed and simulated 
500 mb circulation patterns during four selected 5-day periods when blockings occurred. 
Unlike the major winter blockings over the Pacific and Atlantic, which are the mani- 
festations of zonal wave number n = 1 or 2 (see Kung etal. 1990), the summer blockings 
are more localized and have a smaller scale. This makes identification of summer blocking 
in the zonal wave-number domain difficult. Despite the smallness of the spatial scale, 
however, the periods of the blocking episodes are not necessarily shorter than those of 
winter blocking (see Kung and Baker 1986; Kung etal. 1990). As listed in Table 1, and 
shown in Fig. 2, summer blockings in the Eurasian and Pacific regions are persistent and 
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Figure 3. Observed and simulated 500 mb circulation during four 5-day periods in July and August 1979 
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continually reinforced in the same general regions with only a slight shift in location. For 
instance, in the eastern Pacific from 23 July to 15 August, despite a small movement of 
the block and some weakening and revival of the system, there is, essentially, one 
persistent system during the later period of the simulation. 

4. Forecast skill 

The forecast skill of the 46-day summer simulations for the experiment and control 
are shown in Fig. 4 for the north Pacific sector (30-86°N, 120°E-120°W) and the north 
Atlantic sector (30-86°N, 120-5°W) in terms of anomaly correlations of the 500 mb 

Anomaly Correlation of 500 mb 2 

1.0 
.80 
.60 
.40 
.20 
O 

-.20 
-.40 
-.60 
-.80 
- 1.0 

1 July-15 August 197® 




1 July-15 August 1979 

Figure 4. Anomaly correlations of the 500 mb geopotential-height field in the north Pacific and Atlantic 

sectors during the simulation. 
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eeoDOtential-height field. The sectors are identical to those used in a preceding study 
(Kung er «/ 1990) for the winter 30-day simulations. Both the experiment and control 
indicate a considerable predictability of the 500 mb pattern during early and late July in 
the Atlantic sector, and in mid July in the Pacific sector, when the blocking patterns are 
realistically simulated in respective sectors. For the last 10 days of simuhtion in ear y 
August the experiment with the updated s.s.t. fields shows a very high forecast skill m 
^e Pacific sector, whereas the forecast skill is lost in the control with climatology s.s_ t.s. 
It is noted in Figs. 2 and 3 that the control fails to produce the persistent Pacific blocking 
during this period. As shown in Fig. 1, from late July to August there are changes to the 
s s t anomaly patterns, whereas the patterns show very little change in July. The 
experiment updated with observed s.s.t.s during the simulation is apparently able to 
maintain a high forecast skill in the later stage of the 46-day simulation when the realistic 
eastern Pacific blocking is generated. However, the control run without the benefit of 
updatTng with observed s.s.t.s collapses, as shown by the 10-14 August control patterns 
IFie 31 8 It is noteworthy that, when a realistic blocking is simulated, a high forecast skill 
is obtained in the sector where the blockings develop, and this fact is consistent with the 

winter simulations as reported previously (Kung et al. 1990). . . nt .. 

As the experiment and control show a striking difference in the last period of the 
simulation it will be worthwhile to examine the forecast skill of this period for smaller 
p-’n irp s illustrates the 500 mb anomaly correlation in four longitudinal sectors 
of^the nort he r n hemisphere (3(M20°E, 90°E-180°, 180-90°W and 90-40°W) from 30 to 


Anomaly Correlation of 500 mb Z 
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82°N during the last 15 days of the 46-day simulation. The high forecast skill of the 
experiment is most clearly observed in sectors 180-90°W and 90°E-180°. For the 180- 
90° W sector, we may note that the persistent summer blocking in the eastern Pacific and 
western North America is realistically simulated in the experiment. Some skill is also 
apparent in the 90-40°W sector when North America is under the influence of the Pacific 
blocking. In reference to the high forecast skill in the sector of 90°E-180°, it is noted 
that the experiment also simulated the Eurasian blocking during this period. 

From Figs. 4 and 5 it is seen that the high forecast skill is associated with realistic 
simulation of blocking. This is also the case for the winter simulation in which the high 
skill in the northern hemisphere is due to the high forecast skill of planetary waves 
n = 1 and 2 (Kung et al. 1990). In the case of the summer simulation, however, the 
blocking episodes are localized and the association of blocking simulation and forecast 
skill should be sought in a smaller-scale range than that of the major winter blocking. 
Figure 6 compares the error kinetic energy of long waves n = 1-7 in simulations and 
persistence. The error kinetic energy is defined as the mean-square error of the simulated 
500 mb wind from the observed wind or, in the case of persistence errors, by the mean- 
square difference between the initial and observed wind. For n = 1—4 the error kinetic 
energy of simulations is generally larger than that of persistence, implying that there is 
no appreciable forecast skill of long waves in the summer simulation. However, the error 
kinetic energy of the experiment and control, particularly that of the experiment, becomes 
smaller than that of persistence for n = 5, 6 and 7 but becomes large again for shorter 
wave numbers to be comparable with that of persistence. It appears that the scale range 
of summer blocking occurrences contribute to the forecast skill of the summer simulations. 

Since the appearance of summer blocking is localized, the small error kinetic energy 
of the mid and short waves in the simulation should be related to the blocked area. The 
local error kinetic energy in six longitudinal strips (30-120°E, 120°E-180°, 180-90°W, 
90-40°W, 40°W-0° and 0-30°E) in the latitudinal band of 30-82°N is shown in Fig. 7. 
Examination of Fig. 7, with respect to the times and locations of simulated blockings 
(Figs. 2 and 3), indicates that the small error kinetic energy (i.e. the high forecast skill) 
of the simulations is found in the regions of confluence in the downstream of realistically 
simulated blockings. The relatively low forecast skill in the blocked area is caused by a 
slight shift of the simulated blocking in areas of low error kinetic energy, although a 
realistic simulation of blocking has been achieved in an overall sense. By combining 
longitudes of the downstream area and the blocked area, a realistic simulation of summer 
blocking, such as the August eastern Pacific blocking, implies a high forecast skill over 
nearly half of the northern hemisphere in the mid and high latitudes. 


5. Concluding remarks 

The summer simulation experiment using the high-resolution GLA GCM with 
realistic ocean surface temperatures was able to generate persistent blocking events in 
Eurasia and the eastern Pacific during the 46-day period from 1 July to 15 August 1979. 
The generation of summer blocking episodes appears, at least in part, to be due to the 
heating at the surface boundary over both the oceanic and continental areas. The control 
run with the s.s.t. climatology failed to produce such blockings in the later portion of 
the simulation. It is shown that there is considerable forecast skill in the area of the 
blocking circulations, particularly in the confluence area downstream. It appears that, 
with the updated s.s.t.s, the simulation of realistic summer blocking episodes in the 
northern hemisphere seems possible, at least for a one-and-a-half month period, with a 
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Error Kinetic Energy at 500mb (30° - 82 N) 



ie^ 21 ie 
1 July - 16 August 1979 


for wave number n = 1-7 during the simulation period for 


Figure 6. The 500 mb ^experimental runs. 
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high forecast skill when the blocking pattern develops or is reinforced. Simulation beyond 
this period has not yet been attempted. 

We assume that the skillful predictions of blocking episodes during July in both the 
experiment and control may have been helped by the improved initial data, implying that 
initial conditions include proper energy sources available. Nevertheless, this experiment 
involved only one initial state, and the results should not be generalized. As shown in 
Palmer’s (1988) study, the variability in prediction skill is strongly related to fluctuations 
in the Pacific/North American mode of frequency variability. The generality of the 
reported simulation and the energetics of processes involved in blocking simulation are 
being investigated in a continuation of the study, and will be reported in future articles. 
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ABSTRACT 

It has been established that 
blocking episodes often dominate the 
hemispherical winter circulation and 
local summer circulation of the 
northern hemisphere. A high forecast 
skill of the general circulation 
model (GCM) is recognized with the 
development of blocking. The major 
winter blockings are formed through 
constructive interference of 
traveling wave n-1 and stationary 
wave n-2. The former is controlled 
by barotropic wave -wave interaction 
and the latter by jji situ warming 
over the Pacific and Atlantic. A 
comprehensive examination of extended 
range forecasting with the Goddard 
Laboratory for Atmospheres GCMs is 
offered with 31-day simulations of 
winter blockings and 46-day 
simulations of summer blocking. The 
importance of realistic sea surface 
temperatures (SSTs) in blocking 
simulation is revealed with the daily 
updating of observed SSTs during the 
GCM simulations. 

The empirical long-range 
forecasting scheme at the range of 1 
month to 1 year is formulated with 


multiple regressions for various 
northern hemisphere phenomena and 
monthly- seasonal values of 
temperature and precipitation. A 
critical examination of the 
forecasting experiments is offered. 
The forecasting of Mei-Yu in Taiwan 
indicates that large-scale global 
teleconnections can be used for long- 
range forecasting of regional 
variables with certain modifications. 

1 . INTRODUCTION? 

It is becoming increasingly 
clear that the development of 
numerical forecasting at the range of 
1 month to 1 year requires the major 
task of constructing coupled ocean- 
atmosphere models. For this purpose, 
basic studies are imperative to 
obtain pertinent information on the 
predictability of the atmospheric 
general circulation model (GCM) and 
the role of ocean surface heating in 
model forecasting. Before the 
eventual goal of the long-range 
forecasting capability of the GCM is 
attained, however, the development of 
an empirical-regression scheme for 
practical regional forecasting is 
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also desirable. Actually, studies in 
both areas complement each other for 
the better physical understanding of 
the behavior of the atmosphere. The 
purpose of this paper is to highlight 
our recent and present effort in 
these two areas. 

In the northern hemisphere, the 
winter circulation is often dominated 
by a sequence of blocking episodes 
associated with the activities of 
low-frequency planetary waves. The 
extreme anomalies of summer 
precipitation are often identified 
with summer blocking patterns . For 
the practical importance of the 
phenomena, as well as the known 
extended predictability associated 
with blocking circulation (e.g., 
Shukla 1981; Holloway and West 1984), 
we focused our GCM-related study on 
simulation of blocking. Since the 
generally preferred locations of the 
northern hemisphere winter blocking 
(i.e., the North Pacific and 
Atlantic) are decided by the land-sea 
contrast and topography, it may be 
postulated that interannual 
variations of blocking patterns are 
mainly controlled by sea surface 
temperatures (SSTs) through heat 
release in the ocean-atmosphere 
interaction. Further, the response 
of the GCM to SST anomalies in 
blocking simulation may generate 
pertinent information for the 
development of the joint ocean- 
atmosphere model. In Section 2 of 
this article, the analysis of winter 
blocking circulation is presented in 
the zonal wavenumber domain in 
reference to prevailing SST 
anomalies. Sections 3 and 4 discuss 
the long-range simulation of winter 
and summer blocking episodes with the 
GCM in reference to initial data and 
SST anomaly field. 

The empirical-regression 


approach on the basis ■ of 
teleconnection is presented in 
Section 5. The quantities 

experimented in forecasting include 
the Indian summer monsoon, and 
monthly temperature and precipitation 
in the middle latitude locations. 
Currently a forecasting experiment of 
the Southeastern Asian monsoon (Mei- 
Yu) , In Taiwan, is in progress. 
Special problems associated with 
regional scale long-range 
forecasting, such as for Taiwan Mei-- 
Yu, are presented for examination. 

9 RACKING CIR CULATION IN THE 
WAVENUMBER DOMAIii 

Early in the study of the 
general circulation, Saltzman (1959) 
proposed that large-scale quasi- 
stationary flow systems are 
maintained by a non-linear barotropic 
transfer of kinetic energy from 
smaller cyclone- scale disturbances 
which have baroclinic energy sources. 
More recently, diagnostic studies 
(e.g. Murakami and Tomatsu 1965, 
Paulin 1970) identified blocking 
phenomena with large-scale baroclinic 
conversion. However, Hansen and 
Sutera (1984), by contrasting the 
blocking and non-blocking periods , 
reported the significance of 
nonlinear wave -wave interactions to 
support the kinetic energy of 
blockings. Egger et al . (1986) 

pointed out that the traveling 
weather system interacts with the 
block. Hansen and Chen (1982), in 
studying a case of Atlantic blocking, 
found that this block was forced by 
the nonlinear interaction of Intense 
baroclinic synoptic- scale waves with 
barotropic ultralong waves. 

In our early study (Kung and 
Baker 1986), the spectral energetics 
of the northern hemisphere 
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circulation during winter of the 
First GARP (Global Atmospheric 

Research Program) Global Experiment 
(FGGE) year was investigated with 

gridded analyses of observational 
data and parallel simulation 
experiments. The gridded analyses 
and simulation experiments utilized 
NASA Goddard Laboratory for 

Atmospheres (GLA) 4° latitude x 
5°longitude coarse resolution model 
(see Kalnay et al. 1983) with 
climatological SSTs. The study has 
shown that the characteristic 

energetics features associated with 
the blocking situations in the 
northern hemisphere can be described 
in terms of the barotropic nonlinear 
wave-wave interaction L(n) and 
baroclinic conversion C(n) where n is 
the zonal wavenumber. The pronounced 
single winter blocking is developed 
and maintained by L(l) from the 
kinetic energy source in n-3-10. In 
the case of the double blocking, both 
L(l) and L(2) support the development 
of blocking with L(2) as the largest 
source. The source of kinetic energy 
at n-3-10 is maintained by the 
abundant baroclinic conversion C(n) 
from the available potential energy 
in this spectral range. 

A general description of energy 
transfer in the general circulation 
for blocking development is reported 
by Tanaka and Rung (1988) on three- 
dimensional normal mode functions. 
Figure 1 shows the temporal 
variations of barotropic (jt?- 0 where m 
is the vertical index) and baroclinic 
(m-3-10) energies in the zonal mean 
motion (a?-0) , ultralong waves (n-1-2) 
and synoptic waves (n-3-15) over the 
northern hemisphere during the FGGE 
winter. The time variation of the 
zonal baroclinic energy (n-0 , / 7 ?~ 3 - 10 ) 
indicates clear energy peaks at 16 
and 28 December and 9 January 



Fig. 1 Time change of barotropic 
(m-1) and baroclinic (m-3-10) 
energies for zonal mean motion 
(n-0), ultralong waves (n-1-2), and 
synoptic waves ( n -3-15) over the 
northern hemisphere from 10 
December 1978 through 31 January 
1979. Appearances of typical Rex 
blockings are marked by arrows over 
the time axis . From Tanaka and 
Rung (1988). 


superimposed on the seasonal trend. 
There are subsequent increases of 
baroclinic energy of synoptic waves 
(n-3-15, ro-3-10) in the time series 
through a process of baroclinic 
instability. The time lag is about 5 
days. The time variation of 
synoptic - scale barotropic energy 
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(n-3-15, m-0) is not shown in the 

figure but it is almost in phase with 
the barocl inic energy. The 

barotropic energy of ultralong waves 
(n-1-2, ro-0) increases three days 

later through the wave -wave 
interaction of kinetic energy, the 
process required for amplification of 
ultralong waves to form the blocking. 
The peaks of (n-1-2, m-0) are 

identified with the most intensified 
pattern of blockings, as marked in 
the figure. After the full 

development of the blocking (i.e., 
after reaching the mature stage of 
the blocking) , the baroclinic energy 
of ultralong waves peaks. The time 
series in Fig. 1 further reveals 
evidence of the upscale energy 
cascade from synoptic-scale waves to 
planetary waves. Also, as shown in 
Kung and Baker (1986), under the 
developed blocking situation the 
meridional heat transport increases 
dramatically, and this apparently 
leads to the peak of the baroclinic 
energy of ultralong waves in the 
mature blocking stage. 

Confirming that the blocking 
can be linked to both baroclinic and 
barotropic types of wave activities, 
the monthly circulation patterns of 
the northern hemisphere winter were 
analyzed with daily observations of 
the 500 mb field and monthly sea 
surface temperature (SST) analyses 
during 34 winter seasons from 1955 to 
1989 by DaCamara et al . (1992). The 
blocking index I at longitude X is 
evaluated by 

I(k) = Z(\,62°N) -ZU,42°A0 

where Z is the 500 mb geopotential 
height. This index was originally 
introduced by Lejenas and 0kland 
(1983) and used in Kung et al . (1989, 
1990) with modifications. Positive 


values of I indicate the existence of 
a quasi-meridional dipole which is 
formed by a high-pressure cell 
poleward and a low-pressure area 
equatorward. The quasi -meridional 
dipole is a characteristic shared by 
both diffluent (Rex- type, see Hex 
1950) and meridional (fl-shaped) 
blockings. Additional scans are also 
made for cases of tilted orientation 
of meridional dipoles and cases of 
blocking partially out of the chosen 
latitudinal band. Pacific and 
Atlantic are the two regions of 
preferred winter blocking activities 
in the northern hemisphere. The 
single dominant blockings in the 
Pacific and Atlantic are identified 
by PAC and ATL. The concurrent 
occurrences of blocking in the 
Pacific and Atlantic during the same 
general period constitute the double 
blocking, and are identified by DBL. 

Examination of SST anomalies in 
three subperiods , of approximately 
ten years, reveals that the anomaly 
patterns of these subperiods are 
quite different and can be associated 
with dominant flow patterns of the 
subperiods. The SST anomalies during 
the subperiods are shown in Fig. 2. 
In the subperiod 1955-1966, ATL and 
DBL show high frequencies of 
occurrences, and the Pacific and 
Atlantic both show generally warm 
SSTs . The SST pattern of this 
subperiod may approximate a composite 
of ATL and DBL types of SST anomalies 
strongly leaning toward DBL type due 
to large SST anomalies during DBL. 
The subperiod 1967-1978, with high 
frequency of ATL and NBL, shows a 
general cooling pattern with a 
moderate intensity, and the pattern 
should satisfy the SST conditions for 
ATL and NBL types of flow. The PAC 
flow has the highest frequency during 
the subperiod 1979-1989. The strong 
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without noticeable blocking 
activities. All other months are 
under the dominance or influence of 
one of these three types of blocking. 
Each episode is a dominant and often 
recurring circulation pattern during 
the period. The lengths of the 
periods identified for these episodes 
are from 25 days to 45 days. Figure 
3 shows the time - longitude plots of 
the daily blocking index for two 
episodes in each category. It 
appears that the double "blocking 
events are more intense and 
persistent than single blockings. 
However, the Pacific and Atlantic 
components of double blocking events 
are not exactly in phase with certain 
time lags between maxima in the 
Pacific and Atlantic. This may 
suggest involvement of traveling 
planetary waves. It is interesting 
to note that the persistent blocking 
episodes involve the recurrence of 
blocking patterns in the same general 
locations during the period of 
episodes, strongly suggesting the 
existence of prevailing boundary 
conditions in terms of SST anomalies 
that favor the development of 
blocking . 

As the wave -wave interaction of 
kinetic energy is associated with the 
amplification of planetary waves, it 
is pertinent to explore the 
constructive interference among 
planetary waves for the development 
of a major winter blocking. The 
interference of n-1 and 2 is 
confirmed with all episodes selected 
in this study. Figure 4 

demonstrates, for one episode in each 
of three blocking categories, the 500 
mb trough- ridge diagrams for n-1 and 
2 in the 54°-70°N band. It is 
noteworthy that the primary locations 
of major blocking episodes during the 
course of blocking development are 


the longitudinal segments where n-1 
and 2 interfere constructively. The 
positions of n-2 maximum amplitude 
stay approximately the same in the 
general area of 0° and 180° 
throughout the winter, although the 
strength of the wave is subject to 
time variation. On the contrary, the 
phase angle and amplitude of n-1 are 
both variable, showing that n-1 
travels slowly around the northern 
hemisphere. When the traveling n-1 
constructively interferes with the 
stationary n-2, the major blocking 
develops. If n-1 is relatively weak 
and n— 2 is strong, n— 2 is responsible 
for creating the double blocking 
pattern. This is the case of DBL 
shown in Fig. 4. It is also 
noteworthy that weakening and 
redeveloping blocking in the same 
blocking episodes correspond well 
with the periodic weakening and 
amplifying of n-1 and 2. 

Anomalies of 500 mb 2 and SST 
are averaged for four episodes in 
each blocking category, and shown in 
Fig. 5 for the area between 20° and 
60 °N. Because the northern boundary 
of the area in the figure is 60 °N, 
the positive anomalies of Z in the 
higher latitudes are only partially 
shown for blocking episodes. Yet the 
sharp contrast of positive Z 
anomalies in the north and negative 
anomalies in the south are well 
recognized in the Pacific for PAC, in 
the Atlantic for ATL, and in both 
Pacific and Atlantic for DBL. In 
Fig. 5, it is evident that PAC 
episodes are associated with the 
generally cool SST anomalies in the 
middle - latitude Pacific and western 
Atlantic. For ATL episodes, the 
negative anomalies of SST are found 
in the middle - latitude Atlantic and 
eastern Pacific. However, the 

negative anomalies found in the 
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average pattern for ATL are not as 
distinct as the negative anomalies 
for PAC. Although SST anomalies of 
individual ATL episodes are as 
distinct as those of PAC episodes 
the location of the Atlantic cooling 
is more variable than for PAC 
episodes, making the average Pattern 
less distinct. Comparing the PAC and 
ATL episodes, it is evident that a 
general cooling in the Pacific an 
Atlantic favors the development of 
major single blocking in the Pacific^ 
When the cooling is limbed tro 
Atlantic and part of the Pacifi , 
with a warm sector in the Pacific, a 
single blocking is likely to develop 

in the Atlantic. 

The SST anomalies of DBL in 
Fie. 5 show definite warm SST 
patterns in most areas of the Pacific 

and Atlantic. This sharply contrasts 
with SST anomaly patterns in PAC and 
ATL situations, where the cool SSTs 
are associated with dominant single 
blocking. It is noted with Fig. 4 in 
the preceding discussion that the 
double blocking situation is 
developed and maintained by 
stationary n-2 . The dominant n-2 
should be supported by the stationary 
baroclinic energy source and we 
confirm that there is no barotropic 
energy available for n-2 through 
nonlinear wave -wave interaction. Th 
warm SSTs associated with DBL thus 
indicate that the double blockings 
are developed with in situ warming of 
the northern hemisphere cir ^ u ^ 10n . 
The cold SSTs associated with PAC and 
ATL then can be interpreted to favor 
the nonlinear wave -wave interaction 
to maintain the traveling n-1 that 
couples existing n-2 to form the 
single blocking pattern. Therefore, 
the distinct energetics difference in 
single and double blocking patterns 
emerges. The double blocking pattern 


is baroclinic in nature whereas the 
single major blocking depends on the 
barotropic nonlinear process o 

energy input. Because of the long 
characteristic time of SST variation 
and consequent persistence o 

blocking pattern, the study of 
blocking in relation to SST anomalies 
is relevant in the long-range 
forecasting of the northern 
hemisphere winter circulation. 


3 T.QNG-RANGE simulation of, 

UTNTER BLOCKING 


In our earlier study by Kung 
and Baker (1986) , the simulation of 
blocking by coarse resolution GLA GCM 
with a climatology SST only produced 
a weak Pacific blocking for January 
1979 The failure of the model to 
amplify n-1 to produce a pronounced 
blocking is attributed to the lack of 
L(l) in the simulated circulation. 
The simulated winter circulation 
shows a strong baroclinic conversion 
at the whole spectral range, 
particularly in planetary waves n-1, 

2 and 3. However, these are only 

associated with the formation of a 

weak blocking, and the eddy kinetic 
energy in the planetary-scale range 
cascades down to the short-wave range 
through wave-wave interaction. 

In the following study by Kung 
et al. (1989), four numerical 
simulations of the global atmosphere 
for January 1979 were analyzed to 
study the energetics of blocking 
formation. Both coarse and high 
resolution (2‘x2.5‘) GLA GCMs were 
employed with GLA and Geophysical 
Fluid Dynamics Laboratory (GFD ) 
initial datasets. All simulations 
continued to use the climatology SST 
for January. Among four simulations 
chat were examined, those by high 
resolution GCM tended to generate 
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realistically strong blocking with 
compatible energetics, as in the 
observed blocking episode. It is 
also noted that the simulation with 
the GLA initial data was able to 
amplify n-1 to produce the second 
blocking in the Atlantic in late 
January following the first one in 
the Pacific. The simulations by 
coarse resolution were failures. 
Apparently, the coarse resolution 
model resulted in a downscale energy 
cascade preventing a proper upscale 
input at n-1 . 

As a logical next step, the 
January 1979 global simulations were 
repeated by Kung et al . (1990), with 
the SST anomaly field updated during 
the model integration. Both the 
coarse and high resolution GLA GCM 
are employed with the GLA initial 
dataset. A simulation discussed 


above (Kung et al . 1989), which was 
performed with the high resolution 
GCM utilizing GLA initial conditions, 
is used as the control run to examine 
the SST updating simulations. The 
SST field is provided by a blended 
analysis of in situ (ship and buoy) 
and satellite-retrieved data. The 
SST data retrieved from satellite 
radiances were obtained by the 
multichannel technique of Susskind 
and Reuter (1985) using the high 
resolution infrared sounder and the 
microwave sounding unit on Ti^os-N. 
The approach at GLA is fundamentally 
different from the current 
operational approach at the National 
Environmental Satellite data and 
Information Service in that the 
surface and atmospheric conditions 
are determined from the radiative 
transfer equations with an iterative 


TABLE I 


Three simulation experiments and the control run of the January 1979 elobal 
atmosphere with GLA GCM, and blocking episodes identified in the observed and 
simulated Northern Hemisphere circulation (after Kung et al. 1990); 



Run 

GCM grid 

(latitude x longitude) 

1/1/79 0000 Blocking period 
GMT (day /mo/year) 

initial data 

- 

Observation 



3/1-14/1/79 





12/1-29/1/79 

— 

Exp A 

4°x5 ,> 

GLA 

3/1-16/1/79 

27/1-30/1/79 


Exp B 

2“x2 . 5° 

GLA 

3/1-16/1/79 

18/1-31/1/79 


Exp C 

2°x2 .5" 

GFDL 

7/1-14/1/79 

16/1-25/1/79 

— 

Control 

2°x2 . 5° 

GLA 

3/1-12/1/79 

21/1-25/1/79 

28/1-30/1/79 
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From Kung. DaCamara, Baker, Susskind and Park (1990). 
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scheme. The basic retrieval system 
and analysis technique are detailed 
in Susskind et al. (1984) and 
Susskind and Reuter (1985). The 
monthly mean field of the in situ and 
satellite retrieved data are blended 
to the 4 °x5 ° grid by solving 
Poisson's equation, subject to 
appropriate internal and external 
boundary conditions. This technique 
is similar to the operational blend 
analysis used at the National 
Meteorological Center (NMC) (Reynolds 
1988). 

For the GCM experiments , daily 
gridded SST fields were obtained by 
linear interpolation of the monthly 
mean blended fields, assuming that 
the observed monthly means are 
located in the middle of each month. 
The assumption is valid due to the 
smallness of time variation in the 
SSTs during winter months. Three 
simulation experiments involved in 
this series of study and the control 
run are listed in Table I with the 
identified blocking period. It is 
noted that the high resolution 
version of the GCM used in 
Experiments B and C is identical to 
the coarse resolution version in 
respect to physical 
parameterizat ions . 

Comparing the simulated 500 mb 
circulations with the observed 
circulation during a 5-day period of 
the developed Pacific blocking (first 
blocking) , Exp B shows the closest 
blocking pattern to the observation, 
both in the shape and the location of 
the block. The general pattern of 
the entire northern hemisphere 
circulation also is much better in 
Exps B and C than in Exp A. During a 
similar 5-day period of the second 
blocking in the Atlantic (Fig. 6), 
Exp B again shows a blocking pattern 
in closest agreement to the 


observation, although the development 
of the Atlantic blocking is delayed 
for several days in the simulation. 
Exp C produced a dipole structure. 
However, the blocking is dislocated 
to the east, the jet is more zonal, 
and the blocking period is shorter 
than in Exp B and the observation 
(see Table I). Exp A, using the 
coarse resolution GCM, indicates the 
formation of the Atlantic blocking, 
but the occurrence is substantially 
delayed and the amplitude and 
longitudinal extent of the blocking 
are also limited. 

Examination of simulation blocking 
episodes apparently indicates that 
the simulation by a high resolution 
GCM, using a realistic SST field as a 
boundary condition, is capable of 
generating a second major blocking 
event. As clarified in Kung et al. 
(1989), the high resolution GCM Is 
able to realistically produce the 
first (Pacific) blocking episode 
through the upscale kinetic energy 
input into the ultralong waves by the 
wave-wave interaction. Since the 
development of blocking immediately 
follows the initialization of the 
GCM, we may assume the existence of 
an adequate amount of spatially 
distributed baroclinic energy (see 
Tanaka and Kung 1988) within the 
initial data. For the second 

(Atlantic) blocking event toward the 
end of the simulation period, 
however, the required source of 
baroclinic energy must be provided by 
some other means. This may suggest 
that the realistic ocean surface 
heating has provided an adequate 
baroclinic energy source during the 
simulation, producing a realistic 
second blocking even after the 
effects of initialization are 
presumably lost. 

To examine the forecast skill 
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of the simulations, the anomaly 
correlations of the 500 mb 
eeopotential field in the North 
Atlantic sector (30°-86°N, 100°-5 W) 
are presented in Fig. 7 for Exps A 
and B and the control run* The 
anomaly correlation of Exp B can be 
seen to represent the group of high 
resolution simulations with Exps B 
and C. The figures show that 
following the initial state, there is 
high forecast skill for all three 
simulations in the Atlantic sector. 
Although the skill scores drop in 
mid -January, toward the end of 
January the relatively high skill of 
Exp B becomes apparent with the 
development of a realistic Atlantic 
blocking. Exp A shows some 

insignificant skill at the end of 


January associated with a rather 
limited blocking with the coarse 
resolution GCM. Throughout the 

simulation period, Exp B exhibits a 
relatively high anomaly correlation, 
while the control run, without a 
realistic SST field, shows a marked 
deterioration toward the end of the 
month . 

The forecast skill of ultralong 
waves may be examined by separating 
zonal wavenumbers n-1 and 2 from 
other wavenumbers in the error growth 
of geopotential height and error 
kinetic energy during the simulation 
period. The former is expressed in 
terms of a root-mean square error in 
the 500 mb geopotential height field 
and the latter in terms of a mean- 
square error in the 500 mb wind 


t 

5 0 
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Fie 7 Anomaly correlations of the 500 mb geopotential field in the North 

Atlantic sector during the January 1979 simulation experiments. From Rung, 
DaCamara, Baker, Susskind and Park (1990). 
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field. Both are taken in the 38° - 
70°N latitudinal band, and it is 
indicated that the ultralong waves 
(n— 1-2) in all of the simulations 
show high skill relative to the 
persistence up to 17-18 days. Beyond 
that, Exps B and C, using the high 
resolution GCM, still show good skill 
toward the end of the month. 
However, for shorter waves of n-3-10, 
the forecast skill is only 
recognizable for the initial 7-10 
days of the simulation. Clearly the 
improved blocking simulation with the 
realistic SST field, particularly the 
generation of the second blocking 
event in the Atlantic toward the end 
of the one -month period, is for the 
improved prediction of the ultralong 
waves . 

Although the improved blocking 
simulation has resulted from the use 
of a realistic heating field, we 
cannot determine in this study if it 
is due to the improved initial field 
of SSTs , their realistic time 
variations during the simulation, or 
their combined effect. This may be 
studied in an additional series of 
simulations with various initial 
times. It is important to note that 
only one initial state (1 January 
1979) is involved in this series of 
simulations. As such, the results 
can not be generalized. If tne 
results of future blocking 
simulations indicate the importance 
of updating SSTs during the model 
integration, then one possible 
conclusion would be that a coupled 
ocean-atmosphere model is critical 
for significantly improving extended- 
range forecast skill. Another view 
would be that in these simulations, 
the upper limit of what would be 
possible with a coupled ocean- 
atmosphere model is tested, since 
updating the observed SSTs during the 


atmospheric model integration may 
imply a correct ocean model SST 
prediction. 

4. LONG-RANGE SIMULATION OF 
SUMMER BLOCKING 

Following the preceding 
simulation of winter blocking during 
January 1979, experiments have been 
performed by Kung et al . (1992) for a 
46 -day, summer period from 1 July to 
15 August 1979, to simulate summer 
blocking. The summer simulation also 
utilized the same high resolution GLA 
GCM, During the simulation, SST 
anomalies were updated daily with 
observations, A parallel control run 
used climatology SSTs in lieu of 
observed SSTs. 

The SST anomalies in Fig, 8 and 
the persistent blocking patterns of 
the eastern Pacific in the later half 
of the simulation period in Fig. 9 
seem to indicate that the summer 
blocking in the eastern Pacific is 
also associated with the positive SST 
anomalies in the eastern Pacific. 
The warming persisted and magnified 
during this period, while the 
blocking is strengthened. As the 
scale of summer blocking is smaller 
than that of winter blocking, it is 
consistent with the finding of 
DaCamara et al . (1992) that the 

smaller blocking tends to be more 
baroclinic in nature. This in situ 
association of warm SST anomaly patch 
and persistent summer blocking may 
well suggest that the realistic 
simulation of the August eastern 
Pacific blocking in the experiment is 
for the updating of SSTs during the 
model integration. The control run 
with the climatology SST failed to 
generate this blocking for the lack 
of realistic surface heating. It is 
also shown that there is a 
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considerable forecast skill in the 
area of summer blocking circulation, 
particularly in the confluence area 
of the downstream. Simulation beyond 
this period has not yet been 
attempted. 

s empir km zMSSESgloS forecasting 

Regression forecasting on the 
basis of teleconnection is important 
since it does not base itself on the 
extrapolation of time series, which 
assumes some steady state or steady 
trend of variations. It is an 
application of synoptic analysis in 
the long-range frame. It implicitly 
relates the physical processes, known 
or not yet known, to the future state 
of the atmospheric circulation beyond 


the deterministic range of the 
current GCM forecasting. Thus, as 
stated by Hastenrath (1985) the 
success of regression forecasting 
depends on the successful 
identification of effective 
predictors, which is difficult due to 
physical processes not yet well 
known. It also should be noted that 
the operational mode , or the pattern 
of the general circulation, undergoes 
considerable variation in a decade or 
two, making the established single 
regression useless after a certain 
period. Pointing to these 

difficulties, the logical strategy 
will be the use of the multiple 
regression approach, with predictors 
selected objectively from a large 
group of meteorological and oceanic 


• Recorded 



Fig. 10 Recorded and forecasted dates of the Indian summer monsoon onset 

on the Kerala coast. From 'Kung, Min, Susskind and Park (1982). 
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parameters . 

Figure 10 demonstrates the 
recorded and forecasted dates of the 
Indian summer monsoon on the Kerala 
coast, by Kung and Sharif (1980, 
1982). In their 1980 study, 8 
predictors were utilized, but they 
are all from the April mean upper air 
parameters over the forecast point. 
In the 1982 study, only 5 predictors 
were utilized, but they were 
objectively selected from January to 
April mean upper air parameters over 
India and Australia, SSTs of the 
Indian region, and precipitation in 
Eurasia (representing the snowcover) . 
A marked improvement in the 1982 
study compared to the 1980 study is 
obvious. The forecast gaps of 1966, 
1969, 1972 and 1977 in the 1980 study 
were resolved in the 1982 study. The 
contrast of these two sets of 
forecasts from 1958 to 1977 clearly 
indicates the importance of 
identifying useful predictors, also 
demonstrating the validity of our 
objective method to select predictors 
from the large pool of meteorological 
and oceanic variables. 

Similar experiments were 
performed by Kung and Tanaka (1985) 
to forecast monthly temperature and 
precipitation with antecedent upper 
air parameters and SST with a 
modified scheme of multiple 
regression. Parameters in both the 
northern and southern hemispheres 
were employed. The scheme uses five 
predictors with the database for the 
period 1963-1983. The scheme chooses 
the first predictor on the basis of 
cross-correlation between the 
predictand and variables for the 
general regression: 

V = a 0 + 3 1 X 1 + a 2 X 2 + 

where a is the regression 


coefficient, x is the predictor and n 
is the number of predictors, and n - 
5 in this case. The selection of the 
predictors and fitting of the 
regression coefficients are performed 
simultaneously in this analysis. In 
order to select the first predictor 
x x from m possible predictors, the 
correlation coefficient between y and 
x, r(y, x*), is computed in a time 
series with each of m possible 
predictors. The . first predictor Xj 
is then selected by 

Max[r 2 {y.Xj ) , i 

i.e., the correlation coefficient of 
x x gives the maximum ^(y^x^) among 
the computed m correlation 
coefficients. After selection of x 1 , 
a simple regression is formed by the 
least squares fitting, 

y = a ± x x + c x 

where c l is the first residual 
component and r(x L , cj) - 0. 

The second predictor x 2 is 
selected by examining the correlation 
between the residual c x and the 
possible predictors remaining: 

Max[r- (e 1< x i ) , i = 1, 2, -.m-l] . 

The first residual, e l( is then 
regressed linearly on x 2 as 

e i = W C 2 

where c 2 is the second residual and 
r ( x 2 ,£ 2 )~0- a 2 determined by the 
least squares method again. After 
the selection of predictors x x 
through x n , the regression 
coefficients a 0 through a n , are 
recomputed with a conventional method 
of the multiple regression. 

This scheme is formulated so 
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July Ts In Northeastern U.S. 

(Forecast with February Data) 


• Recorded 
QFortcastfrd 



Correlation between forecasted and recorded surface 
temperature 



Fig. 11 ForecasC experiment of July surface temperature in the 

northeastern United States. From Kung and Tanaka (1985). 
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that the regression residuals are 
minimized with the fewest possible 
predictors. The collinearity among 
predictors is effectively reduced 
even though the correlation among 
selected meteorological parameters 
may not be zero. It is also noted 
that the successive regression with 
residuals is applied only for the 
selection of predictors. Since the 
regression coefficients are computed 
after the selection is completed, the 
present scheme is a multiple 
regression scheme in the conventional 
sense. The comparison of the 
regression coefficients with those 
obtained during the selection, shows 
that their difference is negligibly 
small. This indicates the validity 
of the present scheme. In Kung and 
Tanaka (1985) the scheme gives a lead 
time of predictors a two to eleven 
month period preceding predictands, 
minimizing the collinearity among 
predictors. 

An example of high forecast 
skill in this series of study is 
shown in Fig. 11, forecasting July 
temperatures in the northeastern 
United States. The correlation 
coefficient between the recorded and 
forecasted temperature is 0.92, and 
the mean square error 0.3°C. It is 
also noted in the lower half of Fig. 
11 that the correlation between the 
forecasted and recorded values is not 
necessarily the highest immediately 
preceding the event. In this series 
of experiments, as illustrated in 
Fig. 12, it is suggested that the 
forecasting on the basis of a simple 
predictor is inadequate, and that the 
use of five predictors seems 
sufficient. In the experiments, the 
predictors are selected independently 
every year without involving the 
forecast year's data. However, in 
most cases, the same predictors are 


R 2 MSE 



Number of Predictors (n) 

Fig. 12 Coefficient of multiple 
determination and mean square error 
as functions of number of predictors 
in predicting January surface 
temperature in the northeastern 
United States. From Kung and Tanaka 
( 1985 ). 


selected from the large pool of 
variables. 

Through our model and empirical 
forecast experiments, the importance 
of SSTs as predictor is well 
recognized by the better forecasts 
which result when the SSTs are 
involved. The release of latent heat 
from the ocean surface is the major 
source of energy input to the 
atmosphere, and the variations in 
location and amount of energy input 
essentially determine the deviations 
of the general circulation from the 
average pattern. Many important 
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works exist in this area (e.g. 
Bierknes 1969; Kawamura 1986; Namias 
1959, 1976; Rovntree 1972). We have 
pursued the problem with specific 
interest in establishing large-scale 
teleconnections between predictands 
and SST predictors . 

In Park and Kung (1988), we 
examined the first two principal 
components of the North American 
summer temperature field, and have 
shown that they have considerable 
prognostic value. It is noteworthy 
that the summer temperature of the 
midwestern to eastern part of North 
America is solely determined by the 
first component. Through ^ the 

teleconnection of the principal 
components and preceding SST field, 
it appears that the thermal fields 
over the North Pacific in the 
preceding winter and spring are 
closely related to the North American 
temperature in the following summer. 
As shown in Fig. 13, the negative 
anomalous SST departures in the 
central North Pacific precede the 
higher summer temperature in the 
Midwest, and the positive departures 
precede the lower summer temperature. 
The signal of anomalous SST 
departures appears in the fall and 
reaches its maximum in the winter and 
spring. As this is an interesting 

relationship previously uncovered, it 
also should be pointed out that the 
relationship may offer a convenient 
way to obtain a quick outlook of the 
coming seasons . 

The El N i n o - S o u t h e r n 
Oscillation (ENSO) is known as a 
global phenomenon whose influence is 
often detected in quasi-stationary 
waves as well as in the related 
synoptic- scale weather events in the 
extratropical latitudes during the 
northern hemisphere winter. Although 
Chen (1982) indicated that the forced 


Rossby mode is weak or absent during 
the northern hemisphere summer, the 
signal of ENSO was detected in the 
North American summer temperature in 
our study, suggesting that the 
influence of El Nino on the general 
circulation persists beyond the 
winter half of the year. It is 
unlikely, however, that El Nino has a 
direct influence on the atmospheric 
circulation in the northern 
hemisphere summer at times when the 
anomalous tropical ocean surface 
condition is modified. Instead of 
direct influence, the ocean- 
atmosphere interaction, such as in 
the Aleutian region and the central 
North Pacific in the transition 
seasons, is closely related to El 
Nino, possibly playing a role in 
determining the downstream 
temperature pattern over North 
America in the subsequent summer. 

Our present multiple regression 
scheme, as described in this section 
and detailed by Kung and Tanaka 
(1985) , is also being applied to the 
forecasting of localized phenomena. 
For example, Kinoshita (1989) 
successfully predicted the recess of 
Southeastern Asian monsoon (Mei-Yu or 
Bai-u) with the scheme. At present, 
a collaborative investigation is in 
progress with the Central Weather 
Bureau, Taiwan, R.O.C. to develop 
forecasting schemes for Mei-Yu and 
seasonal rainfall in Taiwan. Since 
our preceding studies of empirical 
forecasting, the database has been 
expanded to include upper air and SST 
data from 1955 to present, and we may 
experiment forecasting beyond the 
defined data period. For instance, 
if we formulate regression equations 
with a 20-year database, we may 
forecast from 1975 to present without 
involving the forecast year during 
the data period. The forecast 
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experiment is thus a completely 
independent experiment in reference 
to basic regression analysis. 

Our ongoing study shows that 
the regression scheme may be applied 
in the restricted-area phenomena in 
Taiwan so long as the phenomena are 
part of large-scale systems, which is 
the case for Mei-Yu and seasonal 


rainfall. However, it also has 
become clear that there are 
considerable noises .of local origin 
affecting the regression forecast. 
The noises seemingly result, at least 
in part, from the interaction of 
large-scale and local-scale 
(subsynoptic and meso scale) effects 
and it is imperative that the 


TABLE II 

Suamary results of 1990 1 ^ 0 !” Forecast..! values represent 

fretted single -regress ton forecasts at various points. Dates of the onset 
and recess are the day from April 1. 


Forecast 


Predictor 


Oct 1989 


Nov 1989 


Dec 1989 


Jan 1990 


Feb 1990 


Onset date 


Recess date Period (days) Rainfall (mm)_ 


700 mb T 
500 mb Z 
SST 

700 mb T 
500 mb Z 
SST 

700 mb T 
500 mb Z 
SST 

700 mb T 
500 mb Z 
SST 

700 mb T 
500 mb Z 


46 . 5 

48.0 

49.5 

51.0 

49.0 

48.5 

45.3 

46.3 

49.7 

47.5 

46.8 

i 

50.0 

51.0 


77.6 

79.7 
77.3 

76.7 

79.0 

76.0 


(a) 


78.7 

81.8 


85.0 


(b) 


(a) 

7 

_(b) 

_(a) 

(a) 


31.5 

34.5 

31.2 

30.3 

30.0 
28.8 

34.4 

34.0 

27.0 

36.8 

37.5 

i 

29.5 
30.3 


480 

447 

454 

452 

499 

485 

437 

356 

471 

494 

481 


(b) 


(b) 


385 

483 


Notes : 


(a) 

(b) 

(c) 


SlI^l'rec^of'lOTo'Tlaim'Si-Vu at Taipei are: 

Onset 48 
Recess 79 
Period 32 days 
Rainfall 420 mm 
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forecasting scheme be properly 
adjusted to suit the local 
situations. Our study indicates 
that, except for years of abnormal 
Mei-Yu (i.e. very early or late 
onset, drought, etc.) the set of 
single regression may be utilized 
effectively if many data points are 
used independently and the cross - 
correlations are properly updated 
every year. The real-time forecast 
of 1990 normal year was thus 
successfully carried out, but the 
abnormal 1991 Mei-Yu (late onset, 
short period) needed a multiple 
regression of three predictors. 
Table II shows the results of 1990 
real time forecast with antecedent 
real-time data from October 1989 to 
February 1990; each value in the 


table representing the average o,f 
single regression forecasts at points 
where strong teleconnections are 
observed between predictands and 
monthly values of predictors. The 
official records of Mei-Yu at Taipei 
are . the onset May 18 (Day 48 from 
April 1), recess June 18 (Day 19) t 
period 32 days, and total rainfall 
420 mm. Thus the forecasts are very 
close to the actual observation. In 
the table it is noted that these 
independent real-time forecasts with 
successive monthly data are very 
consistent and prediction is possible 
from fall. 

Table III compares the 1990 and 
1991 real-time forecasts for Mei-Yu 
onset with single regression and 3- 
term multiple regression. The 




Month of 
real-time data 
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results of single regression 
forecasting are disturbed in 1991, as 
affirmed by the large standard 
deviations in 1991 among independent 
regression forecasts. However, the 
standard deviations are only slightly 
larger for multiple regression in 
1991 than in 1990 when multiple 
regression is used. It may suggest 
the general applicability of the 
multiple regression scheme in both 
normal and abnormal years. It is 
likely that in the abnormal year, the 
operational mode of the general 
circulation deviates from the normal 
prevailing pattern, making the single 
regression less useful. However, 
when the multiple regression is 
utilized, the possible change of the 
operational mode is already 
incorporated in the selection of 
multiple predictors. 

Figure 14 illustrates the 
results of forecast experiments for 
Mei-Yu period and rainfall from 1968 
to 1988 with January predictors. 
Observational data in January of the 
forecast year are used as the real- 
time data. The 3-term multiple 
regressions are estimated with a 10- 
year database, and used for three 
successive years. Thus each year's 
forecasts are done with three 
independently obtained sets of 
regression, each set representing a 
group of multiple regressions at 
different grid points of strong 
cross-correlations. The best value 
among three sets of regressions for 
each year is plotted in the figure. 
The levels of departure among 
forecasts are indicated in Table III. 
The figure clearly demonstrates the 
utility of multiple regressions. 

6 , REMARKS 

Our recent efforts in extended- 


range numerical model forecasting and 
empirical long-range forecasting have 
been reviewed in this paper. The 
model forecasting research is focused 
on the blocking simulation in the 
northern hemisphere winter and 
summer. The blocking appears to 
possess considerable predictability 
and thus useful forecast skill in the 
model integration for at least one 
and one-half months. At present, the 
study continues with specific 
interest in the response of 
simulation to updated SST fields 
during the model integration. The 
study is intended to obtain essential 
information regarding the feasibility 
of constructing the joint ocean- 
atmosphere model. 

The study in empirical 
forecasting for the range of 1 month 
to 1 year is progressing in the 
regression forecast of seasonal 
rainfall in Taiwan. The study is 
intended to provide viable 
information to apply the 
hemispherical teleconnections to 
localized predictands. We intend to 
formulate a generally applicable 
procedure for operational regression 
forecasting in a general setting in 
the northern hemisphere . 
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abstract 


From a 34-year data period, twelve major winter blocking episodes in 
the northern hemisphere are selected for examination: four episodes for each 
‘ h j "“^'. pacific ^Atlantic and double blocking. It is shown that the single 
blockings in the Pacific or Atlantic are formed through constructive inter- 

ference^f the traveling zonal wave of n=I nnd stationer, wave of n=Z The 
ference or me s through the nonlinear wave-wave in- 

n=l is supported by the energy mpm g and Atlantic. T he 

and Atlantic are formed by the 
SnTrJ n=2 when the traveling n=l is weak. The negative anorn^y fields 
of 4a surface temperatures in the Pacific and Atlantic are associated with 
° . . hinrkina in the Pacific and Atlantic. On the contrary, the posi- 

tive anomalies in both the Pacific and Atlantic are associated with concurrent 
Pacific and Atlantic blockings, indicating the baroclinic nature of the double 

blocking. 


1. INTRODUCTION 


The winter circulation in the northern hemisphere is often dominated by a sequence 
j •_ .I,. Pacific and Atlantic regions in association with the activities 
of blocking episoties m tiie d g forecast skill exists for 

of low-frequency planetary waves, u » notea u Holloway and West, 1984; 

^ch P c J s ‘^ cn J^ C ^- m “^ J u Jy of winter blocking particularly important. A major 
Kung et aL, \ 8 ^ transition of the extratropical circulation from states 

^ nCem .^^ t^ilt septic scale disturbances to the low-frequency regimes dominated 
byTasi-stationary planetary waves, and the maintenance of an established blocking pattern. 
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Despite the known difficulty in precisely defining the blocking phenomena, the general 
characteristics of quasi-stationary ridges of blocking circulation are well recognized, and a 
wealth of literature exists on the subject since the early studies of Elliot and Smith (1949), 
Berggren et aL (1949), Rex (1950) and others. More recent development in this area was 
reviewed in Saltzman et aL (1986). During a blocking period, distinct obstruction of the 
northern hemisphere westerlies is apparent in and around the block. However, from the 
viewpoint of the major blocking, it is not a local system. As shown by many in recent years 
(e.g., Hansen and Chen, 1982; Hansen and Sutera, 1984; Kung and Baker, 1986; Kung et 
aL, 1990; Saltzman et al., 1986; Tanaka, et aL 1986; Vakalyuk, 1985), the diagnosis in the 
zonal spectral domain is an effective approach to reveal essential mechanisms involved in 
growth and maintenance of blocking patterns. 

There have been diverging opinions on the energetics nature of blocking among re- 
searchers. Saltzman (1959), in his early study of the general circulation, proposed that the 
large-scale quasi-stationary systems are maintained by a nonlinear barotropic transfer of ki- 
netic energy from smaller cyclonic-scale disturbances which have baroclinic energy sources. 
Diagnostic studies by Hansen and Chen (1982), Murakami and Tomatsu (1965), Paulin (1970) 
and others identified blocking phenomena with large-scale baroclinic conversion. However, 
Hansen and Sutera (1984), by contrasting the blocking and nonblocking periods, reported the 
significant wave-wave interactions to support the kinetic energy of blocking. From the vortic- 
ity consideration, Green (1977), Illari (1984), Shutts (1983) and others also demonstrated that 
the synoptic-scale eddy forcing supports the formation of blocking. Our preceding studies 
of observed and simulated blocking episodes (Kung and Baker, 1986; Kung et aL, 1989) 
revealed that the development of northern hemisphere winter blocking is associated with the 
nonlinear wave-wave transfer of kinetic energy from synoptic-scale disturbances to planetary 
waves. For observed blocking episodes during December 1978 and January 1979, Tanaka and 
Kung’s (1988) energetics diagnosis in three-dimensional normal mode expansion depicted an 
ordinary transfer of energy from the zonal baroclinic component to the barotropic component 
of planetary waves via synoptic scale conversion. Kung et aL (1990) further demonstrated 
that a high resolution general circulation model is capable of producing two successive re- 
alistic blockings in the Pacific and Atlantic through one winter month in January 1979 if 
the sea surface temperature (SST) field is updated daily with observations during the model 
integration. 

As blockings are a group of quasi-stationary circulation patterns of similar synoptic 
appearance, but likely of different causes and mechanisms involved, it is reasonable to find 
diverging diagnostic results among various case studies. To approach a problem of complexity 
such as blocking, both the case study and long-term climatological analysis are necessary. 
In our preceding work (DaCamara et aL, 1991) monthly circulation patterns of the northern 
hemisphere were analyzed with daily 500 mb field and monthly sea surface temperature (SST) 
anomaly field during 34 winter seasons from 1955 to 1989. All 102 monthly circulations 
were classified into 5 categories of patterns according to dominant blocking activities of 
the month. Pacific blocking, Atlantic blocking, double blocking, sequential blocking and no 
blocking. It appeared that, on the gross average, the Pacific, Atlantic and double blocking 
dominated winter circulation and can be associated with characteristic SST anomaly field in 
the Pacific and Atlantic. However, by taking the gross averages of the long-term period, 
variability among different cases can be easily masked. Thus a parallel case study of the 
same time period is highly desirable. 

In this study, prominent major blocking episodes are the focus of our attention. Twelve 
major winter blocking episodes are selected from the daily analyses of 500 mb flow from 1956 
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to 1989: four episodes for each category of Pacific, Atlantic and double blocking. Through 
examination of prominent cases, this study attempts to distinguish the characteristics of these 
blockings in their development and maintenance. Manifestations of the blocking flow pattern 
are examined in terms of the geopotential height and associated anomaly fields. Activities 
of planetary waves are examined by trough-ridge diagrams, kinetic energy and its wave- 
wave interaction to provide energy input at zonal wavenumber n=l and 2. The patterns 
of associated SST anomaly fields are then considered in terms of energy source for in situ 
baroclinic conversion and wave-wave interaction. It is noted that this is a case study and not 
an ensemble analysis, although twelve cases are more than the number of cases usually treated 
in most case studies. We are not attempting a composite climatological analysis of numerous 
blocking episodes during a long-term period. Such an attempt would be difficult, since various 
causes and thus various mechanisms may be involved in many cases of developing blocking. 
Rather, we would like to examine characteristics of major episodes in three categories of 
blocking development which show common prominent appearances in each category. 

2. DATASETS AND SCHEME OF ANALYSIS 

Datasets utilized in this study include the daily northern hemisphere geopotential height 
at 500 mb (Z) in December, January and February from December 1955 to February 1989, and 
monthly SST analyses for the same period. The Z fields are obtained from the daily National 
Meteorological Center (NMQ octagonal grid analyses at 1500 GMT from 1955 to 1957, and 
1200 GMT for the rest of the period. The octagonal grid analyses are edited and bilinearly 
interpolated to a 4° x5° latitude-longitude spherical grid from 18°N to 90°N. The wind fields 
are computed using the geostrophic approximation, and their Fourier transforms are obtained 
at latitude circles 4° apart Monthly SST analyses are acquired from the Comprehensive 
Ocean Atmosphere Data Set (Slutz et al, 1985) for the period 1956-1979, and NMC real- 
time analyses (Reynolds. 1988) for the period 1970-1989. Consistency of both analyses is 
verified for the overlapping period, and a complete time series of 1956-1989 is generated at 
the 2°x2° grid. The Z and SST anomalies at each grid point are computed in reference to 
the 34-year mean value of the month. 

A blocking is recognized when the blocking index I at longitude A is greater than 50m 
along a longitudinal sector: 

/(A) = Z(A,62°iV) - Z(A,46°JV) 

This index was originally introduced by Lejenns and 0kland (1983) and used in Kung et 
aL (1989,1990) with minor modifications. Positive values of I indicate the existence of a 
quasi-meridional dipole which is formed by a high-pressure cell poleward and a low-pressure 
area equatorward. The quasi-meridional dipole is a characteristic shared by both diffluent 
(Rextype, see Rex, 1950) and meridional (fl-shaped) blockings. In the additional scan for 
cases of tilted orientation of meridional dipoles and cases of blocking partially out of the 
chosen latitudinal band, Z(A±5°. 62°N) is substituted for Z(A, 62°N) in computing the 
blocking index I. The blocking identified by I is verified by manual inspection of the flow 
pattern. Pacific and Atlantic are the two regions of preferred winter blocking activities in the 
northern hemisphere. The single dominant blockings in the Pacific and Atlantic are identified 
by PAC and ATL in this paper. The recurrence of blocking in the Pacific and Atlantic during 
the same general period constitute the double blocking, and are identified by DBL. 
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r 1116 kuiet [c energy at wavenumber n, K(n), and the nonlinear transfer of kinetic energy 
from wavenumber m to wavenumber n, L(n,m). are evaluated in reference to planetary waS 
activities. Evaluation of K(n) is described in Saltzman (1957) and our previous studS (e g 
Kung and Baker 1986). The term LOun) for an open system, defined by Two itu HSg* 

domett SU ? aces ’ cx P hc,ll y incorporates the fluxes across the boundary of the 

domain. Its evaluation is discussed by Kanamitsu (1980), Kanamitsu et al (1972)Tand Tsay 

Tl IQQU 9 ^ 8 '^! C . vak3Uon ,n 11118 stU(J y Allows their formulation (see DaCamara It 

damat^on mh S .h! ? m i . (1991) : Sincc our ^ysis * based on the single level 

data at 500 mb, the vertical flux term is not involved in computation. 

3. BLOCKING EPISODES AND PATTERNS OF CIRCULATION 

,, TwelVC . maj0r episodes are selecled from the daily analyses of 500 mb flow for 

T*' ^ f0m epiS0des ’ “<* fa lh '« ‘Rories: PAC. ATL mi 

DBL. These episodes are identified in this paper as PI, P2, P3 P4 for PAC- Al a? ai a a 

for ATL ; and Dl, D2. D3, D4 for DBL. Table 1 lists perij^ 

dcSCnptl0n f of development for each of twelve episodes. Of 108 winter months 
of the 34-year period, we document only 26 months without noticeable blocking activities 
(DaCamara et a/.. 1991). All other months are under foe dominance or ^uenc^of one 
°i bockm 8- Episodes chosen in this study show prominent patterns of 
circulation. However, they are chosen for their prominence and persistence to study 
their distinguished characteristics with least interference of secondary features. 

As described in Table 1, each episode is a dominant and often recurring circulation 
pattern during the penod The lengths of periods identified for these episodes are from 25 
ays to 45 days. Figure 1 shows the time-longitude plots of foe daily blocking index for all 

If hw?‘ SOdeS u ^ Study ‘ Initiation * development, decaying and re-enforcement 

of blocking are identified by manual examination of such time-longitude diagrams. Some 
subjecuve judgement is inevitable in identification, but as shown in Figure 1, identification is 
blocking occurrence is generally recognized when foe contour lines 

Zl TrI T u ,S ^ 0Wn m F,gurc 1 11131 0,(5 winter blockings in foe Pacific and Atlantic 
are formed over foe ^oceanic area, although they may extend into foe continental area. From 

both Figure 1 and Table 1 it appears that foe double blocking events are more intense and 
persistent than single blockings. However, foe Pacific and Adamic components of double 

^ A ^ n C . V r ent l are not exact, y m P ha8e with certain tbne lags between maxima in foe Pacific 
and Atlantic. This may suggest involvement of traveling planetaiy waves It is interesting 
to note that the persistent blocking episodes involve foe recurrence of blocking pattern! in 
foe same genend locations during the period of episodes. This may suggest foe^xistence of 
prevailing boundary conditions that favor foe development of blocking 5 The SST anomaly 
patterns will be examined in this regard later in this paper, ^ 

Figures 2, 3 and 4 are foe 10-day mean circulation patterns and anomaly fields of Z 
between 20 N and the pole, for each episode of PAC, ATL and DBL when foe blockings 
are fully developed. The prominence of foe Pacific blocking in foe northern hemisphere 
c culation is apparent in Figure 2. As foe blocking dominates in foe Pacific, foe Adantic is 
occupied by basically zonal flow. In foe Pacific, the anomaly field is characterized by a sham 

AthS N n f g3UV ® anomaly vea in foe south and positive area in foe north. During dre 
Atlanuc blocking, the patterns in foe Pacific and Adantic are reversed (see Fieure 3) During 

doub " blocki "S u* characteristic flow paueras and associaeTaSttiy Mds a£ 
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• hnth the Pacific and Atlantic. However, examining Figures 2, 3 and 4. it is evident that 

relSto sW anomaly pattern, this seems reasonable, since the oceanic area of the Pacific 
is much more extensive than that of the Atlantic. 

Table 1. Selected major episodes of PaciSc. Allantic and double blocking. 


Case Period Subperiod 
(day/no/yr) (day /mo) 
and Number 
of days 


PI 


P2 


P3 


22/12/77 22-26/12 

-20/1/78 27/12-15/1 

(30days) 16-20/1 

7/12/80 7-11/12 

-5/1/81 12-16/12 
(30days) 17-31/12 


P4 


27/12/82 
-30/1/83 
( 35days) 


17/12/83 

-20/1/84 

(35days) 


I- 5/1 

27-31/12 

1 - 10/1 

II- 15/1 
16-20/1 

21- 25/1 

26- 30/1 

17-21/12 

22- 26/12 

27- 31/12 

I- 5/1 
6 - 10/1 

II- 15/1 
16-20/1 


Primary 
Location In 
42-62*N 


( 1 ) 

( 2 ) 


180-120"W Developing Rex blocking 

180— 120* W Fully developed Rex blocking 

180-120 *W Decaying Omega blocking 

150*E-150*W Identifiable Rex blocking 
150*E-120*W Decaying ridge 
150*E-120*W Fully developed Rex blocking 

120*E-180* Slowly retrograding from (1) to 

150*E-150*W Inverted Omega blocking 

180-120*W Developing Omega blocking 

150*E-150*W Fully developed Rex blocking 
150*E-150*W Decaying Omega blocking 
120-90 *W Developing Omega blocking 

180-120*W Fully developed Omega blocking 

120-90*W Decaying Omega blocking 

180-120*W Omega blocking 

180— 120*W Fully developed Rex blocking 

180-120*W Decaying ridge 

120-90*U Developing ridge 

180-120*W Developing Omega blocking 

180-120*W Fully developed Omega blocking 

180-120*W Decaying Omega blocking 


( 2 ) 


Al 

26/1- 

26-30/1 


24/2/66 

31/1-4/2 


(30 days) 

5-14/2 

15-19/2 

19-24/2 

A2 

31/1- 

31/1-4/2 


24/2/69 
(25 days) 

5-14/2 

15-19/2 

19-24/2 

A3 

27/12/71 

27/12-10/1 


-20/1/72 

11/1-15/1 


(25 days) 

16-20/1 

A4 

11/1- 

11-15/1 


14/2/79 

16/1-4/2 


(35 days) 

5-14/2 


30*U-30*E Developing ridge 

0* Developing Omega blocking 

6Q*y_0 # Fully developed Rex blocking 

60*W-0* Omega blocking 

60*U-0* Decaying ridge 

60*U— 0* Developing Rex blocking 

6G*V-0* Fully developed Rex blocking 

90— 30*W Fully developed Omega blocking 

90— 30*W Decaying ridge 

E Fully developed Rex blocking 

E Fully developed Rex blocking 

Fully developed Rex blocking 

Developing ridge 
Fully developed Rex blocking 
Slowly retrograding from (1) to (2) 
Fully developed Omega blocking 


30*W-30 < 

30*W-60" 

0-60*E 

60*W-0* 

(1) 60*W-0* 

(2) 90-30*W 
60*W-0* 
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D1 


D2 


D3 


DA 


6/1- 

6-10/1 

(i) 

180-120'W 

4/2/57 


(2) 

30*W-30*E 

(30 days) 

11-20/1 

(1) 

180-120'W 



(2) 

30*W-30*E 


21-30/1 

(1) 

180-120'W 



(2) 

0-60'E 


31/1-4/2 

(1) 

180-120W 



(2) 

0-60'E 

7/12/61 

7-11/12 

(1) 

180-120'W 

-5/1/62 


(2) 

30'W-30'E 

(30 days) 

12-21/12 

(1) 

150'E-150'W 



(2) 

0-60'E 


22-26/12 

(1) 

120'E-180* 



(2) 

90-60 'W 


27-31/12 

(1) 

120'E-180* 



(2) 

60*W-0' 


1-5/1 

(1) 

150'E-150'W 



(2) 

60'W-0 

22/12/62 

22-26/12 

(1) 

180-120'W 

-A/2/63 


(2) 

30'W-30'E 

(A3 days) 

27/12-10/1 

(1) 

150'E-150'W 



(2) 

60'W-0* 


11-25/1 

(1) 

180'-120'W 



(2) 

30*W-30'E 


26-30/1 

(1) 

180-120'W 



(2) 

60'W-0* 


31/1-4/2 

(1) 

150'E-150'W 



(2) 

60“W-0* 

26/1- 

26-30/1 

(1) 

180'-120'W 

29/2/68 


(2) 

30'W-30'E 

(35 days) 

31/1-4/2 

(1) 

120*E-180* 



(2) 

60*W-0* 


5-9/2 

(1) 

120-90'W 



(2) 

60*W-0' 


10-19/2 

(1) 

120-90'W 



(2) 

60'W-0* 


20-29/2 

(1) 

120-90'W 



(2) 

30“W-30'E 


Omega blocking developing in 
both (1) and (2) 

Fully developed Rex blocking 
in both (1) and (2) 

Omega blocking in both 
(1) and (2) 

Decaying Omega blocking in 
both (1) and (2) 

Developing Omega blocking 
in both (1) and (2) 

Fully developed Omega blocking 
in both (1) and (2) 

Fully developed Rex blocking 
in both (1) and (2) 

Omega blocking in both 
(1) and (2) 

Decaying ridges in both 
(1) and (2) 

Omega blocking in both 
(1) and (2) 

Fully developed Rex blocking 
in (1) and (2) 

Omega blocking in (1) and 
Rex blocking in (2) 

Fully developed Omega blocking 
in both (1) and (2) 

Rex blocking in (1) and 
Omega blocking in (2) 

Developing Omega blocking in 
(1) and developing ridge in (2) 
Rex blocking in (1) and 
Omega blocking in (2) 

Omega blocking in both 
(1) and (2) 

Rex blocking in both 
(1) and (2) 

Omega blocking in (1) and 
Rex blocking in (2) 


Under the dominance of a single blocking in the Pacific or Atlantic, the amplification 
of the planetary wave n=l is most obvious in the northern hemisphere circulation (Figures 2 
and 3), whereas the amplification of n=2 occurs in the double blocking circulation (Figure 4). 
This situation may permit the treatment of major winter blockings in the zonal wavenumber 
domain despite their traditional, synoptic treatment as local manifestations. As the wave- 
wave interaction of kinetic energy is associated with the amplification of planetary waves, it 
is pertinent to conjecture that the constructive interference among planetary waves may be 
a factor for the development of a major winter blocking. Austin (1980) reported that n=l 
and 2 tend to interfere constructively in the Atlantic sector, and n=2 and 3 interfere in the 
acific sector. In our preceding study (Kung et al., 1990), it was suggested that during the 
development of a Pacific blocking in January 1979, n=l and 2 came in place at the location 
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Fig. 4. As in Figure 3, but for DBL. 
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of the block. Toward the end of the same month the interference of n=l and 2 in the Atlantic 
was identified with the development of a major Atlantic blocking. 

The interference of n = l and 2 is confirmed with all episodes selected in this study. 
Figure 5 demonstrates for one episode in each of three blocking categories the 500 mb 
trough-ridge diagrams for n=l and 2 in the 54° -70° band. Although the blocking index is 
defined between 46° and 62°N to account for the general pattern of circulation, the 54°-70°N 


P4 



A3 



03 



Fig. 5. Trough-ridge diagrams of n=l and 2 in the 54°-70°N band with 500 mb 
Z during blocking episodes of P4 (PAC), A3 (ATL) and D3 (DBL). The 
contour interval is 50 m, and negative contour lines are dashed. 
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band is used in this application because it is the latitudinal band in which the blocking pattern 
is maximized. It is noteworthy that the primary locations of major blocking episodes (as 
indicated in Table 1) during the course of blocking development are the longitudinal segments 
where n=l and 2 interfere constructively. Madden (1983) found that the interferences of 
stationary and traveling waves of the same longitudinal scale cause time variations in the 
large-scale circulation. It is noted in Figure 5 that the positions of n=2 maximum amplitude 
stay approximately the same in the general area of 0° and 180° throughout the winter, 
although the strength of the wave is apparently subject to time variation. On the contrary, the 
phase angle of n=l, in addition to the amplitude, is variable, indicating that n=l travels slowly 
around the northern hemisphere. It is seen here that, when the traveling n=l constructively 
interferes with the stationary n=2, the major blocking develops. If n=l is relatively weak 
and n=2 is strong, n=2 is responsible for creating the double blocking pattern, as is the case 
of D3 shown in Figure 5. It is also noteworthy that weakening and redeveloping blocking in 
the same blocking episodes, as identified in Table 1, correspond with the periodic weakening 
and amplifying of n=l and 2. 

A comprehensive energetics description is not the purpose of this paper. Neither is it 
possible with the dataset of this study to compute energy transformations, except for nonlinear 
interactions at the 500 mb level. Thus only the kinetic energy K and nonlinear transfer L 
are presented in Figures 6 and 7 to describe characteristic meridional distribution of K(n) for 
n=l,2 and 3 for PAC, ATL, and DBL episodes and associated L(njn) in reference to n=l. 
The characteristics of circulation are first observed by K(n) as kinetic energy and represent 
the strength of circulation. Examination of L(n) is pertinent for blocking cases, since the 
development and maintenance of blocking are closely associated with the nonlinear transfer 
of kinetic energy, which provides the barotropic energy source (Kung and Baker, 1986; Kung 
ef al, 1989). As shown in Figure 6, K(l) of PAC is most pronounced among the three 
categories of blocking, and is also least unstable from episode to episode. For ATL, K(l) 
is less pronounced than for PAC, which is consistent with the observation that the Pacific 
blocking is more dominant than the Atlantic blocking (see Figures 2 and 3). For DBL 
episodes, K(2) is most pronounced among planetary waves, as expected from the dominance 
of the n=2 pattern. 

The L(l), the input of kinetic energy to n=l from all other waves, and L(l,2), the 
input of kinetic energy from n=2 to n=l, are most pronounced for PAC, as seen in Figure 
7. This indicates that the traveling n=l (see Figure 5) has its source at n=2 in forming the 
Pacific blocking. In the case of DBL, n=l is also maintained by n=2 through L(l,2), but the 
magnitude of K(l) is much less than for PAC. For ATL the barotropic energy source of n=l 
is not at n=2. An earlier case study of an ATL episode by Kung and Baker (1986) shows that 
during the development stage of an atlantic blocking n=l is supported by nonlinear interaction 
with synoptic scales smaller than n=3. It is also noted in this study, though not shown in 
Figure 7, for all three categories of blocking situation the wave-wave energy input into n=2 is 
very small: negative for PAC and ATL and only slightly positive for DBL. Thus the energy 
source at n=2 should come from the baroclinic conversion, both for the maintenance of n=2 
itself and in support of transformation to n=l. 

4. SEA SURFACE TEMPERATURE ANOMALIES AND BLOCKING 

From the foregoing discussion of planetary wave activities, it may be questioned if the 
three categories of major winter blocking are associated with different types of energetics 
forcing. It is generally recognized that the SST anomaly fields exercise a major control over 
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Fig. 6. Meridional distributions of kinetic energy K(n) for n=l, 2 and 3 in units 
of m 2 s -2 for four episodes and their mean for PAC, ATL and DBL. 


the winter circulation through heat release in the ocean-atmosphere interaction (e.g., Kung 
et oL 1990; Namias. 1959, 1978; Wallace et aL, 1990). The patterns of SST anomalies are 
examined in reference to the major blocking episodes in this study. 

Anomalies of Z and SST are averaged for four episodes in each blocking category, 
and shown in Figure 8 for the area between 20° and 60°N. Blocking area includes both the 
positive area of Z anomalies in higher latitudes and negative area in lower latitudes. Because 
the northern boundary of the area in the figure is 60°N. the positive anomahes of Z m the 
higher latitudes are only partially shown for blocking episodes (see Figures 2.3, and 14). Yet 
th? sham contrast of positive Z anomalies in the north and negative anomahes in the south 
are well recognized in the Pacific for PAC, in the Atlantic for ATL and in both Pacific and 
Atlantic for DBL. It is noted that the areas of blocking both in the Pacific and Atlantic tend 
to appear in lower latitudes for DBL than for PAC and ATL. 

In Figure 8, it is evident that PAC episodes are associated with the generally negative 
SST anomalies in the middle-latitude Pacific and western Atlantic. For ATL episodes die 
negative anomalies of SST are found in the middle-latitude Atlantic and eastern Pacific. 
However, the negative anomalies found in the average pattern of ATL in Figure 8 are not as 
distinctive as the negative anomalies for PAC. Yet, as shown in Figure 9, SST anomalies of 
individual ATL episodes are as distinct as those of PAC episodes. As the location of the 
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Fig. 7. Meridional distributions of kinetic energy input at n=l by nonlinear wave- 
wave interaction L(l), and its component from n=2 and 3, L(l,2) and 
L(l,3), in units of 10 -5 m 2 s -3 for composites of PAC, ATL and DBL. 


negative Atlantic SST anomalies is more variable than for PAC episodes, the average pattern 
becomes less distinct. This seems to be related with the fact that the circulation pattern of 
ATL episodes is more variable than that of PAC episodes, as evident in the circulation patterns 
of Figures 2 and 3 and latitudinal distributions of K(l) and K(2) in Figure 6. Comparing the 
PAC and ATL episodes, it is apparent that general negative SST anomalies in the Pacific and 
Atlantic favor the development of major single blocking in the Pacific. When the negative 
pattern is limited to the Atlantic and part of the Pacific, with positive anomalies in the Pacific, 
a single blocking is likely to develop in the Atlantic. 

The SST anomalies of DBL in Figures 8 and 9 show definitive positive patterns in most 
areas of the Pacific and Atlantic. This is a sharp contrast with SST anomaly patterns in PAC 
and ATL situations where the negative SST anomalies are associated with dominant single 
blocking. It is noted in the preceding section that the double blocking situation is developed 
and maintained by the stationary n=2 (Figure 5). The dominant n=2 should be supported by 
the baroclinic energy source since no barotropic energy is available forn=2 through nonlinear 
wave-wave interaction. The positive SST anomalies associated with DBL indeed indicate 
that the double blockings are developed with in situ wanning of the northern hemisphere 
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Fig. 8. Anomaly patterns of 500 mb Z and SST for composites of episodes for 
PAC, ATL and DBL. The contour interval is 50 m for Z and 0.5°C for 
SST. Negative contour lines are dashed. 
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circulation. The negative SST anomalies associated with PAC and ATL then can be interpreted 
to favor the nonlinear wave-wave interaction to maintain the traveling n=l that couples 
existing n=2 to form a single blocking pattern. The distinct energetics difference of single 
and double blocking patterns shown by these cases are interesting. The prominent double 
blocking pattern is baroclinic in nature whereas the occurrence of single major blocking 
depends on the barotropic nonlinear process of energy input. 

5. CONCLUDING REMARKS 

Because of the long characteristic time of SST variation and consequent persistence of 
blocking pattern, the study of blocking in relation with SST anomalies may have relevance 
in the long-range forecasting of the northern hemisphere winter circulation. With twelve 
prominent blocking episodes during the 34-year data period, this study has shown that the 
major single blockings in the northern hemisphere winter circulation are developed m the 
Pacific or Atlantic trough constructive interference of traveling n=l and stationary n-2. The 
former is supported by the nonlinear wave-wave interaction of kinetic energy, and the latter 
by the in situ warming over the Pacific and Atlantic. The double blockings situation in 
the Pacific and Atlantic is formed by n=2 when the traveling n=l is weak. These planetary 
wave activities are seemingly associated with the accompanying SST anomaly patterns during 
these major blocking episodes. The negative SST anomalies of the Pacific and Atlantic are 
associated with PAC and ATL, whereas the in situ warming in both the Pacific and Atlantic 

is associated with DBL. .... . , . 

For the purpose of focusing on prominent major blocking situations, the approach taken 

in this study is that of case examination with only twelve major blocking episodes. As such, 
the results discussed only apply to major blocking episodes studied in this paper, and cannot 
be generalized to other blocking circulations. However, the results of this case study may offer 
useful information in considering the formation and maintenance of the northern hemisphere 
winter blockings. At present, a composite diagnosis of all winter blocking episodes during 
the 1956-1991 period is in progress in terms of planetary wave activities and will be reported 

subsequently. 
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1 . 


Introduction 


It is well recognized that the ocean surface provides a major boundary 
forcing to the general circulation of the atmosphere through the process of 
ocean-atmosphere heat transfer. A systematic investigation of the sea surface 
temperature (SST) anomalies in reference to the tropospheric circulation is 
thus a desired step toward better understanding of climatic fluctuations. 

This study is one such attempt by means of principal component analysis to 
identify prevailing modes of variations in SSTs and the tropospheric 
circulation. 

The utility of the principal components, which are also referred to as 
empirical orthogonal functions or eigenvectors, in analyzing the large-scale 
variations of the circulation field has been demonstrated by Kutzbach (1967, 
1970), Kidson (1975a, b), Weare et al . (1976), Weare (1977), Barnett (1978), 
Heddinghaus and Kung (1980), Trenberth and Paolino (1981), Kawamura (1984, 
1986), Park and Kung (1988) and others. One goal of such an analysis is to 
reduce a large number of variables into a manageable set of components, while 
retaining the maximum variance of the original variables, enabling us to 
concisely describe patterns of interrelationship among observed variables. In 
the analysis of highly correlated meteorological fields, a limited number of 
principal components may effectively represent the fundamental modes of 
variations which are highly loaded on these components. It is important to 
note that the empirical orthogonal functions have no predetermined forms, and 
only depend on the interrelationship within the dataset of analysis. Thus, 
they are particularly suitable to investigate anomalous fields of the general 
circulation, for which no known analytical form exists because of complex 
boundary conditions and various scales of nonlinear interactions. 


68 


In this study, pattern vectors and time coefficients are obtained for 
major principal components of monthly mean fields of the global SSTs and the 
Northern Hemisphere geopotential height Z(m) at 700, 500 and 300 mb. 

Variances of intra-annual seasonal variations are eliminated by defining the 
anomaly fields of monthly means in reference to the long-term means of 
individual months for the entire data period from 1955 to 1992. The 
characteristics of the first three components of SST fields are identified and 
discussed. The quality of SST data changes significantly during the data 
period because of the availability of satellite observations in the later 
years. The heterogeneity of data in time series is examined, and its 
elimination is considered. The intra-annual seasonal variations of SSTs are 
also presented with principal components of long-term means of twelve 
individual months. The association of SST components and Z components is 
studied on the basis of inter-annual variations of these components. Cross- 
correlations of SST components with the observed Z fields are further studied 
to examine the seasonal - range predictability of the tropospheric circulation 
in reference to SST anomalies. 
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2 . 


Data and scheme of analysis 


Datasets utilized in this study include monthly global SST analyses for 
the period of 1955-92, and daily northern hemisphere Z analyses at 700, 500 
and 300 mb for the same period. Monthly SST analyses were obtained from two 
different sources : the Comprehensive Ocean-Atmosphere Data Set (COADS) (Slutz 

et al. 1985) for the period of 1955-79 and the NMC SST real-time analyses 
(Reynolds 1988) for the period of 1970-92. Both analyses are on 2°x 2° 
latitude - longitude grids, but the grid points of two analyses are staggered. 
The consistency of both analyses was verified for the overlapping period of 
1970-79, and the COADS analysis of 1958-69 was interpolated on the grid points 
of the Reynolds analysis to make a complete time series of 1955-92. Missing 
data points in early COADS analysis were interpolated in space from the 
surrounding grid points. The SST data were used for the oceanic areas between 
40° S to 60°N. The monthly Z fields were computed from the daily National 
Meteorological Center (NMC) Northern Hemisphere octagonal grid analyses at 
1500 GMT from 1955-57, and at 1200 GMT for the rest of the period. 

The general mathematical procedures for computation of orthogonal 
pattern vectors and their coefficients are as provided in Kutzbach (1967). 

The grid values of monthly mean fields are standardized to zero mean and unit 
variance so that each grid point can have normalized significance in 
describing the spatial variation. As stated in the introduction, the anomaly 
fields are obtained with respect to the multi-annual means of respective 
monthly mean fields for the study of inter-annual variations. This allows for 
the elimination of the seasonal cycle. On the other hand, the intra-annual 
seasonal cycle is presented by subjecting the multi-annual means of 12 
individual months to principal component analysis. The variable matrix, which 
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has normalized departures from time means, is used to compute the 
corresponding correlation matrix, whose elements are the correlation between 
grid points. The eigenvectors of such a correlation matrix represent a set of 
normalized departure fields. The eigenvectors are not rotated in order to 
retain the maximum variance in the first few components (see Walsh and Richman 
1981). The characteristic pattern vector is described as the eigenvector 
multiplied by the standard deviation of its corresponding principal component. 
Thus, characteristic patterns represent the simple correlations between 
respective components and meteorological fields. 

Most of the available principal component analyses of large-scale 
meteorological fields indicate that the first three components are sufficient 
to describe the large-scale variations (e.g., Kutzbach 1970; Weare et al . 
1976; Weare 1977; Park and Rung 1988). As discussed by Overland and 
Preisendorfer (1982), it is important to determine if the principal components 
obtained in the analysis are statistically significant, so that the 
geophysical interpretation of results will be meaningful. As shown in pattern 
vectors of major components (Figs. 2, 4, 5, 6 and 7), large-scale 
characteristic patterns are described mostly by large correlation coefficients 
of 1 0 . 4 1 to | 0 . 8 | , which are significant at the 98 to 99.9% confidence level 
for the time period of datasets. The higher components also show some 
significance in limited areas, but the level of significance in large areas 
drops sharply. In this paper, the discussion is limited to the first three 
components . 
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3. Principal components of SST fields 

The inter-annual variations of coefficients of 1st principal components 
of January and July SSTs for the period of 1955-92, as shown in Fig. 1, 
indicate a large gap at 1970. This phenomenon is also observed in all other 
months. If the SST data for the period of 1970-92 are subjected to principal 
component analysis, the 1st and 2nd components of the new time series match 
almost exactly with the 2nd and 3rd components of the original 1955-92 time 
series. The discontinuity of the 1955-92 data at 1970 in the 1st principal 
component is important as the 1st component contributes 18% of the total 
variance, and is larger than any other component's contribution. This gap 
indicates a major systematic difference of SST data before and after 
(including) 1970, most likely reflecting both the inclusion of satellite data 
and a different scheme of analysis in NMC data from that in earlier COADS 
records . 

In this study, the original 1st component in the 1955-92 dataset is 
disregarded, and the 2nd, 3rd and 4th components are treated as the true 1st, 
2nd and 3rd component of the 1955-92 period. The principal components to be 
discussed hereafter refer to the true components after elimination of the 
original 1st component as a systematic bias. This is justified by the near 
exact match of the 1st and 2nd components of the 1970-92 dataset with the 2nd 
and 3rd components of the 1955-92 datasets. Since we are only interested in 
the first few major components, small-scale variations that are associated 
with the original 1st component, other than the gap at 1970, do not enter into 
the picture and can be ignored. 

Figure 2 shows the pattern vectors for the first three components of 
January and July SSTs, and their corresponding coefficients are shown in 
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Fig. 3. The percentage variances of the first five components out of the 
total 35 components, listed in Table 1, indicate a comparable contribution of 
the first few components with previously available principal component 
analyses (see citations in Section 2). The pattern vectors of the 1st 
components, as shown in Fig. 2, represent the largest-scale spatial variation 
in the Pacific and Atlantic. The variation in the northwest- southeast 
direction is conspicuous in the Pacific with a maximum positive area in the 
tropical region. The 2nd component shows a well defined positive area in the 
eastern tropical Pacific, which is similar to the 1st component previously 
obtained in the regional analyses of the Pacific SSTs (e.g., Weare 1986; Weare 
et al . 1976). Much of the other areas in the Pacific and Atlantic are 
negative, contrasting the large positive areas of the 1st component. The 3rd 
component shows variations in the Pacific and Atlantic in the northwest- 
southeast direction. The pattern vectors beyond the 3rd component show a 
localized, fragmental pattern, and are not included in our discussion of 
large-scale SST variations. 

It is apparent that the pattern vectors of the 2nd component represent 
the El Nino- Southern Oscillation (ENSO) type of SST anomalies. Since the 1st 
component also has a general positive area in the central and eastern tropical 
Pacific, it is interesting to examine the time variations of these two 
components. As marked for January SSTs in Fig. 3, the prominent ENSO winters 
of 1957-58, 1965-66, 1972-73, 1976-77, 1982-83 and 1991-92 are characterized 
by the amplification of both 1st and 2nd components in phase. The combined 
coefficients of the 1st and 2nd components in Fig. 3 are weighted means by the 
respective percentage variances as listed in Table 1. In cases of the 
moderate to weak ENSO winters, either the 1st or 2nd component is amplified, 
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but not both of them in phase. The following summer of an ENSO winter is not 
necessarily associated with the ENSO type SST anomalies. When the 1st and 2nd 
principal components are reversed following the ENSO winter (see Fig. 3), the 
southeastern Pacific is occupied by extensive negative areas. Such is the 
case for July 1964, 1973, and 1988 when negative SST anomalies developed in 
this area. 

When we subject the long-term (1955-92) SST means of the 12 months of 
the year to the principal component analysis, we obtain pattern vectors and 
time variations of coefficients as shown in Fig. 4 for the 1st and 2nd 
components. The variation of pattern vectors are concentrated in the tropical 
latitudes, and the time variations of both the 1st and 2nd components reflect 
the time scale of seasonal variations. Comparing the pattern vectors of 4 
separate months in Figs. 2 and 5, it is seen that basic patterns of 3 
principal components seem to hold throughout the year, but are subject to the 
seasonal changes which are not negligible in terms of strength and local 
modification. The intra-annual seasonal variations presented in Fig. 4 should 
be considered when the persistence of prevailing SST anomalies is to be 
projected. The pattern vectors of January and July, as shown in Fig. 2, are 
modified between winter and summer to exhibit spring and fall patterns, as 
shown in Fig. 5 for April and October. 
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4. Principal components of the Northern Hemisphere circulation 

Pattern vectors of the first three components of the Northern Hemisphere 
circulation in terms of January and July Z fields are shown in Figs. 6, 7 and 
8 respectively at the 700, 500 and 300 mb levels. A similarity in pattern is 
observed through the troposphere from 700 mb to 300 mb, and the 2nd component 
shows the best similarity at all three levels with "centers of action" in the 
same quadrants. These centers of action are stable and dominant, and 
consistent with the generally recognized winter flow patterns. As the first 
two components of winter Z fields determine the general pattern of the winter 
Northern Hemisphere flow, it is also noted in Table 1 that the contribution of 
percentage variance in the winter flow by the 2nd Z component is much greater 
than those by the 2nd components in SST and summer Z. The 3rd component, 
among the first three components, shows the more localized patterns, which is 
a generally observed character of higher components in the principal component 
analysis . 

When the July patterns are compared to the January patterns in Fig. 6-8, 
it is apparent that the patterns of all summer components are weakened at all 
levels, and the scale of dominant patterns becomes generally small. The 
distinguished patterns of the 2nd and 3rd components in winter become 
fragmented in summer with the sharpest contrast at 500 mb. However, 
components in summer show general consistency at different levels. The 
contrast between the summer and winter patterns merely reflects the difference 
between the dominant winter circulation and the significantly weakened summer 
circulation in the Northern Hemisphere. 

Explicit association of physical processes with separate components of 
the circulation fields is beyond the scope of the principal component 
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analysis. However, because of the lower boundary forcing of SSTs, it is 
useful to examine the correlations between components of SSTs and Z fields. 
Table 2 lists the correlation coefficients between the first three components 
of SSTs in January, April, July and October, and those of monthly Z fields at 
700, 500 and 300 mb during the period of 1955-92. The correlation coefficient 
| 0.35 | is significant at the 95% level, and |0.44| at the 99% level. In Table 
2 it is evident that the 1st SST component is strongly correlated with the 1st 
Z component. The correlation is large at 300 mb and 500 mb throughout the 
year, particularly at 300 mb. However, at 700 mb it is limited to the winter 
season. In other words, the 1st components of SSTs with the Z components is 
more significant in the upper troposphere than in the lower troposphere. The 
correlation of the 2nd and 3rd SST components with Z components is generally 
low. Although some sporadic correlations exist in the table, they do not 
indicate an organized pattern of correlation. According to the correlations 
in Table 2, only the 1st SST component shows a significant, persistent 
correlation with the 1st Z component in the mid and upper troposphere. 

As the 1st and 2nd SST components seem to determine the occurrence of 
ENSO events (Fig. 3), it will be interesting to examine if the principal 
components of the upper-air circulation show a specific ENSO mode as in SSTs. 
Figure 9 illustrates the composite 700 mb and 500 mb flow pattern and anomaly 
fields for the prominent ENSO January months of 1958, 1966, 1973, 1977, 1983, 
and 1992. The 700 mb anomaly field resembles the pattern vectors of the 2nd 
principal component of the January 700 mb field (Fig. 6), and the 500 mb 
anomaly field shows a similarity to the pattern vectors of the 3rd principal 
component of the January 500 mb field. It is noted, however, that the inter- 
annual variations of such principal components are not clearly associated with 
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winter ENSO occurrence. The situation is thus different from that of the time 
variations of SST components in association with ENSO occurrence (see Fig. 3). 
Although SST s show a clear large-scale mode of boundary forcing in ENSO 
winters, the responding upper-air circulation does not show such clear ENSO 
components. The anomalous circulation patterns, as often reported for ENSO 
winters, seem to be modifi- 
cations to the basic field of circulation rather than dominant specific large- 

scale patterns . 
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5 . 


Cross-correlation patterns 


In a lag-correlation analysis, Kawamura (1986) demonstrated that time 
coefficients of the first component of the monthly mean North Pacific SST 
tends to lead the 500 mb Z in the extratropical latitudes. The lead time is 
up to 12 months or more, depending on the month of SST. Park and Kung (1988) 
showed that the SST anomalies for preceding seasons in the North Pacific are 
closely related to the first component of surface temperatures in North 
America in the following summer. The latter determines the summer temperature 
in the midwestern and mid-Atlantic United States. It is possible that the 
principal components of SST fields may be useful in the long- (seasonal-) 
range empirical forecast. 

Figures 10, 11 and 12, respectively, illustrate the cross - correlation 
patterns between the 1st, 2nd and 3rd principal components of January SSTs and 
the January, April and July anomaly fields of 500 mb and 300 mb Z. All three 
SST components are well correlated to the 500 mb and 300 mb Z field of the 
same month, January. The contour spacing is 0.2, and the correlation in the 
figure is significant at 90, 95, 99 and 99.5% with the coefficients of |0.28|, 

| 0.33 |, | 0.42 | and |0.45|. In April, with a lead time of 3 months, all three 

SST components show similarities with January flow patterns. Some variations 
from the January patterns, particularly the strengthening of the 2nd component 
correlation in April is noted. The correlations become weaker and more 
fragmented in July, apparently reflecting the weaker and more fragmented 
circulation patterns of the summer. However, the centers of action as 
recognized in earlier months in Pacific, Atlantic, North America and Eurasia 
are still traceable. 
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Comparing cross -correlations of January SST components with 500 mb and 
300 mb flow patterns in Figs. 10-12, it is noted that the characteristic 
correlation patterns, as observed with 500 mb fields, are identified in 300 mb 
Z fields in terms of both spatial distribution and strength of cor-relation 
coefficients. This suggests that the response of tropospheric circulation to 
SST anomaly fields is consistent in the mid and upper troposphere. 

The cross -correlation patterns of the 1st, 2nd and 3rd principal 
components of July SSTs with the July, October and January 500 mb and 300 mb Z 
fields are illustrated in Figs. 13, 14 and 15. The cross -correlation patterns 
observed in July carry through October and January. The patterns observed in 
500 mb are also clearly observed in 300 mb. One clear distinc-tion from the 
January SST cross-correlation is the strong correlation pattern of July SSTs 
with January 500 mb and 300 mb Z fields with a lead time of at least 6 months. 
It is also noteworthy that the cross - correlation with the January Z fields is 
stronger than with the October Z fields. These reflect strong winter 
circulation patterns, and also indicate a useful seasonal - range predictability 
of winter upper-air circulation patterns with the anomaly fields of summer 
SSTs. 


79 



6. Concluding remarks 

As seated in the introduction, the purpose of this paper is to identify 
prevailing modes of variations in SSTs and the tropospheric circulation. The 
results obtained in this study offer useful information in further diagnosis 
of the SST and circulation fields, and also the utilities of the SST anomaly 
field in seasonal-range forecasting of the tropospheric circulation. The 
principal component analysis of monthly mean fields of the global SSTs and the 
Northern Hemisphere tropospheric circulation for the 38-year period as 
presented in this study seems to yield the following specific information. 

1 a change in the source of SST data results in a large systematic 

bias in a long-term SST dataset. However, the bias may be effectively 
corrected if an appropriate principal component is identified and eliminated 
from the dataset. 

2. Periods of strong ENSO correspond to both the 1st and 2nd SST 
components being in phase. 

3. The first few principal components of monthly SST and Z fields are 
sufficient to describe the SST and tropospheric circulation patterns. 

4. Basic patterns of the first three principal components of SSTs 
remain similar throughout the year, but are subject to seasonal changes, which 
are not negligible. 

5. The first SST component is strongly correlated with the first 

component of tropospheric Z field. 

6. Cross -correlations of SST components with the tropospheric 
circulation suggest a useful predictability of the tropospheric circulation in 
reference to SST anomalies. The lead time observed in this study is up to 6 
months for the winter circulation in reference to summer SSTs. 
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Table 1. Percentage variance of the first five principal components 

during the period of 1955-92. 


Component 

SST 

Z(700mb) 

Z(500mb) 

Z(300rab) 

J an . 

July 

Jan . 

July 

Jan . 

July 

Jan. 

July 

i 

17.1 

12.2 

17.5 

17.1 

17.7 

20.0 

20.0 

22.1 

2 

7.4 

8.6 

15.5 

9.9 

13.0 

8.7 

13.2 

12.0 

3 

5.6 

5.5 

9.3 

8.9 

10.3 

8.5 

8.9 

7.0 

4 

5.0 

5.2 

8.7 

7.3 

9.0 

6.8 

7.6 

5.8 

5 

4.4 

4.9 

6.4 

5.9 

6.3 

5.6 

6.7 

5.0 

Total 

39.5 

36.4 

57.4 

49.1 

56.3 

49.6 

56.4 

51.9 
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Table 2. Correlation coefficients between principal components of monthly 
SSTs and principal components of monthly Z fields at 700, 500 and 300 mb during 
the period of 1955-92. 


SST (1st) January April July October 


z 

Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

700 mb 
(1st) 

- .48 

- .63 

- .60 

-.53 

.25 

.44 

-.28 

- .34 

.22 

- .06 

- .29 

.51 

(2nd) 

. 46 

.03 

-.19 

- .06 

- .45 

.24 

-.27 

.02 

.43 

.62 

- .19 

.09 

(3rd) 

- .18 

.11 

- .20 

- . 35 

.40 

- .24 

- .04 

.50 

.13 

.43 

.31 

.05 

500 mb 
(1st) 

-.36 

- .55 

.48 

- .49 

.43 

- .61 

- .45 

- .45 

.24 

- .26 

- .12 

-.51 

(2nd) 

- .43 

- .04 

- .18 

- .23 

- .01 

.24 

- .08 

.15 

- .13 

.08 

- .04 

.12 

(3rd) 

.09 

- .05 

-.23 

- .18 

- .24 

-.17 

.13 

.13 

.06 

-.47 

.56 

.12 

300 mb 
(1st) 

. 56 

.71 

.75 

.65 

.75 

- .84 

.64 

.60 

.61 

.34 

.25 

.49 

(2nd) 

.22 

- . 15 

.16 

.18 

.17 

.14 

- .17 

- .05 

.11 

.13 

- .23 

.04 

(3rd) 

-.21 

- . 16 

.12 

- .18 

.25 

- .01 

.05 

- .19 

.09 

.36 

.42 

- .05 


SST(2nd) 


January 



April 



July 


October 


z 

Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

700 mb 
(1st) - 

.27 

- .46 - 

.19 

.17 

.02 

.03 

- .20 

- .06 

.02 

- .06 

-.16 

.02 

(2nd) 

.20 

- .05 - 

.08 

. 37 

.17 

-.11 

.20 

.06 

- .30 

.20 

- .29 

.18 

(3rd) 
500 mb 

.01 

- .29 

.08 

.15 

- .09 

- .16 

-.02 

- .01 

.17 

.08 

.05 

- .23 

(1st) - 

.06 

- .22 

.04 

- .01 

.01 

- .02 

-.33 

- .34 

.19 

.35 

.21 

.00 

(2nd) - 

.11 

- .10 - 

.07 

.34 

- .07 

.06 

.04 

.09 

.42 

- .04 

-.21 

- .33 

(3rd) 
300 mb 

.45 

.49 

.00 

- .33 

- .15 

- .03 

- .27 

.17 

-.10 

- .18 

-.10 

.07 

(1st) 

.09 

.23 

.12 

-.15 

- .22 

.21 

.25 

.39 

.26 

-.19 

- .11 

.06 

(2nd) - 

.01 

. 11 

.13 

- .25 

.04 

.12 

- .13 

-.14 

- .03 

.02 

- .09 

- . 34 

(3rd) 

.05 

.44 

.11 

- .26 

.25 

.14 

.43 

- .04 

.43 

- .11 

- .01 

.04 

SST(3rd) 


J anuary 



April 



July 


October 


Z Jan 

Feb Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

700 mb 
(1st) - 

.32 

. ii 

.03 

.24 

- . 32 

- .03 

- .03 

-.00 

-.20 

.15 

-.19 

.23 

(2nd) - 

.61 

.02 - 

.01 

- .01 

- .29 

.40 

.43 

.08 

-.04 

- .21 

- .20 

.00 

(3rd) - 

500 mb 

.01 

.04 

,30 

- .35 

.11 

.18 

- .18 

- .34 

-.26 

- .13 

.01 

-.18 

(1st) - 

.26 

.04 

11 

.28 

- .42 

.14 

- .16 

- .18 

.17 

.08 

- .02 

-.19 

(2nd) 

. 64 

.04 - 

04 

- .05 

.05 

- .22 

.01 

- .25 

.12 

- .16 

- .23 

-.15 

(3rd) - 

300 mb 

.03 

- . 17 

18 

- .08 

.25 

-.13 

-.27 

- .17 

.09 

- .03 

.19 

.04 

(1st) - 

.04 

- .22 - 

08 

- . 17 

- .24 

.04 

- .09 

-.22 

- .29 

- .05 

-.01 

.12 

(2nd) - 

.66 

- .04 - 

04 

.12 

.03 

- .05 

.07 

- .28 

-.03 

.11 

.17 

- .23 

(3rd) - 

.01 

- .20 

01 

.01 

- .18 

- .07 

- .04 

- .01 

.17 

.08 

.15 

.17 
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FIGURES 


Fig. 1 Inter-annual variations of coefficients of the principal components 1, 
2 and 3 for the original January and July SST datasets during the 
entire data periods of 1955-92 and the partial data period of 1970-92. 

Fig. 2 Characteristic patterns (pattern vectors) for the first three 

principal components of January and July SSTs during the period of 
1955-92. Positive areas are shaded. The contour interval is 0.2. 

Pig 3 Inter-annual variations of coefficients of the first three principal 
components of January and July SSTs during the period of 1955-92. 

Dots denote years of major winter ENSO episodes. 

Fig. 4 (a) Characteristic patterns for the first two components of long- 

term (1955-92) mean monthly SSTs, and (b) corresponding coefficients 
of principal components. Positive areas are shaded. The contour 
interval is 0 . 2 . 

Fig. 5 Characteristic patterns for the first three principal components of 

April and October SSTs during the period of 1955-92. Positive areas 
are shaded. The contour interval is 0.2. 

Fig. 6 Characteristic patterns of the first three principal components of 

700 mb Z field in the Northern Hemisphere for January and July. The 
contour interval is 0.2 with real lines for zero and positive, and 
broken lines for negative. 

Fig. 7 As in Fig. 6 , but for 500 mb. 

Fig. 8 As in Fig. 6 , but for 300 mb. 

Fig. 9 700 mb and 500 mb January flow patterns and anomalies for composites 

of prominent January ENSO cases . 

Fig. 10 Cross-correlation patterns of coefficients of the first 

January SST component and 500 mb and 300 mb Z fields in January, 

April and July. The contour interval is 0.2 with real lines for 
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zero and positive, and broken lines for negative. 

Fig. 11 As in Fig. 10, but for the second January SST component. 

Fig. 12 As in Fig. 10, but for the third January SST component. 

Fig. 13 Cross - correlation patterns of coefficients of the first July SST 
component and 500 mb and 300 mb Z fields in July, October and 
January . 

Fig. 14 As in Fig. 13, but for the second July SST component. 

Fig. 15 As in Fig. 13, but for the third July SST component. 
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Cross-Correlations of 2nd January SST Component with Upper Air Anomalies 
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Figure 11 


















Cross-Correlations of 3rd January SST Component with Upper Air Anomalies 
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Figure 14 














Cross-Correlations of 3rd July SST Component with Upper Air Anomalies 
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Figure 15 
















